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I  Objectives 

The  research  program  for  FY88  is  characterized  by  two  major  themes,  namely,  i)  the  contin¬ 
uation  of  the  research  effort  in  supersonic  shock  wave-turbulent  boundary  layer  interaction, 
and  il)  the  initiation  of  a  significant  research  activity  in  3-D  hypersonic  shock  wave-turbident 
boundary  layer  interaction.  The  basic  four  objectives,  outlined  in  the  original  proposal  [1] 
and  described  below  are  pertinent  for  both  supersonic  as  well  as  h]rpersonic  research  activity. 

1-1  Develop  and  Validate  Theoretical  Model(s)  for  3-D  Shock  Wave- 
Turbulent  Boundary  Layer  Interaction 

One  of  the  fundamental  steps  in  the  understanding  of  complex  3-D  shock  wave-turbulent 
boundary  layer  interactions  involves  the  development  and  validation  of  accurate  theoret¬ 
ical  models.  Niunerical  effort  has  progressed  from  the  simplified  calculatioiu  of  the  pre¬ 
supercomputer  years  [2,3]  to  the  current  complex  simulations  requiring  significantly  more 
computer  resources.  The  Reynolds-averaged  full  three-dimensional  compressible  Navier- 
Stokes  equations  [4]  are  typically  utilized  with  turbulence  incorporated  through  a  model 
such  as  the  two-equation  Jones-Launder  [5]  or  the  algebrsdc  Baldwin-Lomax  [6]  model  This 
complexity  is  necessary  to  adequately  resolve  the  several  fluid  dynamic  effects  that  play  crit¬ 
ical  roles  in  3-D  turbulent  interactions  including  boundary  layers,  shock  waves,  turbulence 
and  viscous-inviscid  interactions.  Recently  published  results  utilizing  the  above  two  eddy 
viscosity  models  have  shown  overall  good  agreement  with  experimental  observations  on  sim¬ 
plified  3-D  geometries  such  as  the  sharp  fin  [7]  and  the  swept  compression  comer  [8,9].  This 
validation,  achieved  by  comparison  with  experimental  observations,  is  a  continuous  element 
of  the  research. 

The  success  of  the  theoretical  approach  has  led  to  the  inception  of  a  new  activity,  that 
of  modeling  three-dimensional  hypersonic  shock  wave- turbulent  bound2ury  layer  interaction 
(’3-D  hypersonic  turbvdent  interactions’)  in  FY88.  This  represents  an  extension  of  previous 
activity  which  was  focused  on  supersonic  flows. 

1-2  Determine  the  Physical  Structure  of  3-D  Turbulent  Interactions  for 
Selected  Geometries 

Closely  related  to  the  task  of  developing  and  validating  the  theoretical  model  is  that  of 
determining  the  structure  of  such  flow  fields.  Over  the  past  several  years,  significant  progress 
has  been  achieved  in  obtaining  an  overall  picture  of  the  flowfield  for  the  sharp  fin  and  the 
swept  comer.  The  author  has  collaborated  with  C.C.  Horstman  (NASA  Ames  Research], 
B.Shapey  and  S.Bogdonoff  (Princeton  University)  in  this  endeavor.  Based  on  the  analysis 
of  flows  past  the  sharp  fin  [7]  and  the  swept  comer  [9],  the  overall  flowfield  is  observed  to 
be  inviscid  and  rotational  over  a  major  portion  of  the  interaction  domain  except  within  a 
narrow  region  adjacent  to  the  surface  where  the  effects  of  turbulent  mixing  are  significant. 
For  both  configurations  (Figs.  1  and  2),  the  principal  flow  feature  is  a  large  vortical  stmcture 
aligned  with  the  comer  in  agreement  with  the  flowfield  modeb  of  Token  [10]  and  Kubota  and 
StoUery  [11].  A  three- dimeiuional  surface  of  separation  (Fig.  1)  emanates  from  the  line  of 
coalescence  (separation),  and  spirab  into  the  vortical  center.  A  second  surface,  emanating 
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from  upstream,  intersects  the  plate  at  the  line  of  divergence  (attachment),  and  defines  the 
extent  of  the  fluid  entrained  into  the  vortical  structure. 

For  the  current  year,  the  primary  goals  are  (Section  11-2): 

•  To  reaffirm  the  validity  of  the  theoretical  model  consisting  of  the  Reynolds-averaged 
three-dimensional  Navier-Stokes  equation  in  mass  averaged  variables  in  conjunction 
with  the  algebraic  turbulence  model  of  Baldwin  and  Lomax.  This  aspect  is  closely 
linked  to  Research  Task  No.  1  above. 

•  To  provide  three-dimensional  flow  visualization  to  explore  the  interaction  flow  structure. 

•  To  examine  numerically  the  validity  of  the  quasiconical  free  interaction  principle. 

1-3  Investigate  Methods  of  Control  and  Modification  of  3-D  Turbulent 
Interactions 

The  understanding  of  the  flowfield  structure  of  3-D  turbulent  interactions  provides  the  oppor- 
ttmity  for  investigation  of  methods  for  control  and  modiflcation  of  the  interaction  flowfields. 
Specifically,  large  vortical  structures  such  as  those  displayed  in  Figs.  1  and  2  are  clearly 
imdesirable  in  typical  applications  such  as  inlets  where  3-D  turbulent  interactions  are  com¬ 
mon.  Improvements  in  the  control  of  such  interactions  may  lead,  for  example,  to  improved 
aircraft  performance.  Continuing  research  efforts  (Section  11-3)  focus  on  bleed  (suction)  and 
generation  of  opposing  vorticity  by  mounting  another  symmetric  fin  on  the  plate  (Fig.  3). 

1-4  Examination  of  Inviscid  and  Viscous  Effects  in  3-D  Turbulent  Inter¬ 
actions 

Computed  results  have  indicated  that  the  flow  structure  for  the  sharp  fin  and  the  swept  comer 
are  not  only  remarkably  similar  but  also  that  the  calculated  flowfields  are  relatively  insensitive 
to  the  turbulence  model  employed.  This  conclusion  is  derived  from  the  fact  that  two  different 
modeb  successfully  predicted  overall  flow  features  even  though  eddy  viscosity  values  from 
the  two  models  differed  significantly  in  the  major  portion  of  the  boundary  layer  outside  of 
the  viscous  sublayer  repon.  The  fourth  research  task  is  therefore  to  analyse  viscous-inviscid 
effects  in  such  interactions  to  further  enhance  understanding  of  the  inherent  physics. 
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II  Research  Accomplishments  for  FY88 

The  research  program  follows  the  basic  tasks  outlined  in  the  original  proposal  submitted  to 
Or.  James  Wilson  on  4  April  1986.  Four  specific  research  tasks  were  outlined  for  FY88  [12], 
namely,  1}  Development  and  validation  of  theoretical  models  for  3-0  shock  wave-turbulent 
boundary  layer  interactions,  2)  Oetennination  of  physical  structure  of  3-0  turbulent  interac¬ 
tions  for  specific  geometries,  3)  Investigation  of  methods  for  control  and  modification  of  3-0 
turbulent  interactions  and  4)  Examination  of  inviscid  and  viscous  effects  in  3-0  turbulent 
interactions. 

The  overall  research  program  is  closely  coordinated  with  the  computational  research  of 
Or.  C.C.  Horstman  (NASA  Ames  Research  Center)  and  the  experimental  research  of  Prof. 
S.  Bogdonoff  (Princeton  University).  Ore.  Knight,  Horstman  and  Prof.  Bogdonoff  meet 
frequently  to  select  specific  3-0  shock  wave-turbulent  boundary  layer  interactions  for  detailed 
study.  Separate,  independent  computations  are  then  performed  for  the  configuration  by 
Drs.  Knight  and  Horstman,  and  the  calculated  flowfields  are  relatively  compared  with  the 
experimental  measurements  of  Prof.  Bogdonoff.  This  close  collaboration  between  theory  and 
experiment  is  a  critical  element  of  the  success  of  the  program. 

For  the  purposes  of  the  following  discussion,  it  is  useful  to  note  the  general  orientation  of 
the  physical  (z,  y,  z)  axes: 

Streamwise  direction:  z 
X  to  plate:  y 
Spanwise  direction:  z 

II-l  Research  Task  No.  1:  Develop  and  Validate  Theoretical  Models  for 
the  3-D  Shock  Wave-Turbulent  Boundary  Layer  Interactions 

A  new  activity  in  modeling  three  dimensional  hypersonic  shock  wave-turbulent  boundary 
layer  interaction  (’3-0  hypersonic  turbulent  interactions’)  is  initiated  during  FY88.  This 
represents  an  extension  of  previous  activity  which  was  focused  on  supersonic  flows.  A  princi¬ 
pal  objective  of  this  research  is  the  determination  of  the  sensitivity  of  the  computed  flowfields 
in  3-0  hypersonic  turbulent  interactions  to  the  turbulence  model  employed.  Recent  exper¬ 
imental  and  theoretical  results  at  Mach  3  have  shown  that  simple  turbulence  models  (i.e., 
zero  equation  and  two  equation  eddy  viscosity  models)  are  capable  of  accurate  simulation 
of  a  variety  of  3-0  turbulent  interactions  [7,8,9],  while  failing  to  accurately  reproduce  the 
limiting  2-D  turbulent  interaction  [13].  The  success  in  analysing  flows  past  the  sharp  fin 
and  the  swept  comer  is  generally  attributed  to  the  nature  of  the  interaction.  Specifically, 
the  flowfields  investigated  are  observed  to  be  essentially  inviscid  and  rotational,  except  in  a 
small  fraction  of  the  turbulent  boimdary  layer  adjacent  to  the  surface  where  the  effects  of 
turbulence  are  significant.  The  absence  of  reverse  flow  may  also  be  a  factor  [14]. 

A  combined  experimental  and  computational  research  project  has  been  initiated  with  Dr. 
C.C.  Horstman  (NASA  Ames  Research  Center)  to  examine  the  3-0  hypersonic  turbulent  in¬ 
teraction  generated  by  a  sharp  fin.  The  selection  of  this  configuration  is  motivated  by  several 
factors.  First,  through  the  recent  collaborative  experimental  and  theoretical  effort  involving 
the  principal  investigator,  NASA  Ames  Research  Center  and  Princeton  Gas  Dynamics  Lab¬ 
oratory  [9]  a  detailed  understanding  of  the  flowfield  structure  of  the  3-0  sharp  fin  turbulent 


6 


interaction  at  Mach  3  has  evolved.  Second,  a  series  of  3-D  sharp  iin  experiments  are  planned 
at  hypersonic  speeds  (Mach  7  and  10)  at  NASA  Ames  Research  Center  in  FY89^.  Efforts 
to  date  have  focused  on  the  determination  of  the  grid  for  the  calculations.  The  activity  for 
FY89  is  discussed  in  Section  IIL 

II-2  Research  Task  No.  2:  Determination  of  Physical  Structure  of  3-D 
Turbulent  Interactions  for  Specific  Geometries 

As  mentioned  earlier,  the  CTirrent  understanding  of  the  flowfield  for  both  the  sharp  fin  and 
the  swept  comer  is  that  it  is  inviscid  and  rotational  over  a  majority  of  the  interaction  domain, 
except  within  a  narrow  repon  adjacent  to  the  surface  where  the  effects  of  turbulent  mixing 
are  significant,  and  dominated  by  a  single  vortical  structure. 

Kimmel  [15]  performed  a  series  of  experiments  to  investigate  the  effect  of  the  shock 
generating  geometry  on  the  interaction.  Three  configurations  were  selected,  namely,  the 
sharp  fin,  swept  compression  comer  and  the  semicone.  The  geometrical  characteristics  were 
selected  to  obtain  similar  shock  strengths.  Specifically,  the  sharp  fin  angle  a  =  17.5  deg, 
the  angle  of  compression  and  sweep  for  the  swept  comer  a  =  30  deg  and  A  =  60  deg,  and 
the  semiconit  half  angle  7  =  20  deg  generate  equal  strength  shocks  in  the  free  stream.  Each 
configuration  was  examined  at  Reynolds  numbers  of  1.6  x  10^  and  9.1  x  10*^.  For  the  lower 
Reynolds  number,  detailed  surface  pressure  and  kerosene-lampblack  flow  visualization  was 
performed.  At  the  higher  Reynolds  number,  experimental  data  include  both  stirface  pressure 
and  kerosene-lampblack  flow  visualization,  and  boundary  layer  profiles  of  pitot  pressure, 
static  pressure  and  yaw  angle. 

The  three  configurations  exhibited  conical  similarity,  with  the  flow  surface  and  flowfleld 
variables  exhibiting  two  distinct  re^ons.  The  first  region,  extending  from  upstream  to  a 
point  downstream  of  the  shock  system,  displayed  a  ’quasi-conical  free  interaction’.  In  this 
region,  for  example,  the  experimental  pressure  distribution,  displayed  along  a  line  normal  to 
the  inviscid  shock,  was  virtually  identical  for  the  three  configurations  despite  their  significant 
geometrical  differences,  A  second  region,  downstream  of  the  first  region  and  denoted  the 
’model  dominated  region’,  showed  no  similarity  of  profiles. 

The  explanation  for  the  quasi-conical  free  interaction  is  not  clear.  It  is  anticipated  that  all 
three  flowfields  display  a  large  vortical  stmcture.  The  sharp  fin  at  a  =  17,5  deg  is  intermediate 
to  a  series  of  computed  flows  (a  =  10  and  20  deg)  which  have  displayed  a  large  vortical 
stmcture  [7].  The  swept  compression  comer  at  (a.  A)  =  (30, 60)  deg  is  close  to  the  previous 
computation  of  Knight  et  al  [8]  at  (a.  A)  =  (24, 60).  The  semicone,  however,  has  not  been 
computed  previously.  It  is  unclear,  however,  why  this  vortical  stmcture  should  yield  nearly 
identical  pressure  distributions  in  the  quasi-conical  free  interaction  region. 

To  achieve  the  primary  objectives  of  this  task,  both  the  17.5  deg  sharp  fin  and  the  (30, 60) 
swept  comer  have  been  computed  in  FY88.  The  results  are  briefly  described  below.  Further 
details  may  be  fotmd  in  [16]. 

’Hotftman,  Private  communication 
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II-2.1  Results  for  the  17.5  deg  Fin 

The  results  are  first  compared  with  experiment.  Boundary  layer  survey  data  is  available  at 
sixteen  (16)  locations  exhibited  in  Fig.  4.  The  computed  and  experimental  surface  pressure 
are  displayed  in  Fig.  5  at  z,  =  9.6tfoo  with  z«  denoting  the  distance  from  the  model  apex 
along  the  inviscid  shock  trace.  Overall  features  of  the  measured  oistribution  are  predicted 
by  the  computation.  Upstream  influence  (determined  by  initial  pressure  rise)  is  underpre¬ 
dicted  in  the  computation.  The  plateau  region,  observed  by  several  previous  investigators, 
between  separation  and  attachment  is  adequately  reproduced  by  the  computation  though  the 
flatness  of  the  plateau  region  is  not  well  resolved.  The  computed  and  experimental  pressure 
rise  after  the  shock  wave  and  the  downstream  pressure  level  differ  by  less  than  7%.  The 
shallower  computed  pressure  gradient  is  assumed  largely  due  to  inadequate  grid  resolution  in 
the  interaction  region. 

Typical  results  of  the  experimental  and  computational  boundary  layer  profiles  of  pitot 
pressure,  yaw  angle  and  static  pressure  at  station  12  are  shown  in  Fig.  6.  Good  agreement 
is  observed  in  the  boundary  layer  for  pitot  pressure,  yaw  angle  and  static  presstire.  The 
computed  pitot  pressure  does  not  resolve  the  small  overshoot  associated  with  the  compression 
fan  outside  the  boundary  layer.  The  computed  static  pressure  continues  to  increase  toward 
the  inviscid  region  while  the  measured  data  slightly  decreases. 

The  measured  surface  streamline  patterns  for  the  17.5  deg  fin  are  displayed  in  Fig.  7. 
The  computed  result  is  shown  in  Fig.  8.  The  computation  predicts  qualitatively  all  of  the 
main  features  revealed  by  the  experiment.  A  line  of  coalescence  is  formed  outside  of  the 
inception  zone  near  the  model  apex.  The  lines  of  upstream  influence  and  coalescence  are 
approximately  straight  rays  from  a  virtual  origin.  The  computation  imderpredicts  the  angle 
of  the  line  of  upstream  influence  and  line  of  coalescence  by  14.5%  and  15.5%  respectively. 
This  underprediction  has  also  been  observed  previously  [7]  for  other  fin  angles.  A  divergence 
line  is  seen  near  the  fin  which  divides  the  fluid  flowing  over  to  that  being  entrained  into  the 
vortical  structure.  The  surface  streamlines  between  the  line  of  coalescence  and  the  line  of 
divergence  (attachment)  converge  asymptotically  with  the  line  of  coalescence.  No  secondary 
separation  [17]  is  observed. 

The  flow  structure  is  examined  with  exhaustive  sets  of  particle  traces.  For  the  purpose  of 
brevity,  only  a  few  are  selected  for  presentation.  Fig.  9  displays  traces  of  particles  released  at 
the  upstream  end  of  the  computed  domain  at  y  =  0.5^oo-  For  illustration,  the  line  of  coales¬ 
cence  is  also  displayed.  The  streamlines  roll  into  a  large  vortical  structure  with  particles  near 
the  fin  surface  pitching  down  and  away  from  the  fin  beneath  the  vortical  core.  The  orienta¬ 
tion  of  the  vortical  structure  is  coimterclockwise  as  viewed  from  the  base  of  the  fin  leading 
edge.  Particles  released  from  upstream  at  y  =  l.ldgo  are  displayed  in  Fig.  10  with  the  line 
of  coalescence.  The  fluid  particles  at  this  height  are  not  entrained  into  the  vortical  structure 
and  simply  flow  over  the  interaction.  Additional  extensive  particle  trues  are  consistent  with 
these  features. 

II-2.2  Results  for  the  (30,60)  Swept  Compression  Comer 

The  surface  pressure  at  z,  =  lO^oo  of  the  (30,60)  swept  comer  is  displayed  in  Fig.  11.  The 
computed  results  underpredict  upstream  influence.  The  plateau  pressure  level,  the  rapid  pres¬ 
sure  rise  after  the  shock  wave  and  the  sharp  peak  (“comer  pressure”)  near  the  model  comer 
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are  accurately  predicted.  The  slight  drop  at  the  end  of  the  plateau  region  is  not  adequately 
resolved.  The  location  of  the  maximum  pressure  predicted  by  computation  is  upstream  of 
that  measured.  This  indicates  that  the  extent  of  the  computed  results  underpredict  the  size 
of  the  interaction  region.  The  computed  iinal  pressure  level  is  about  10%  below  the  measured 
value. 

The  computed  and  measured  contour  plots  for  the  flowfield  static  pressure  is  shown  in 
Fig.  12.  Overall,  good  agreement  is  obtiuned  between  computation  and  experiment.  Al¬ 
though  the  slope  of  the  upstream  foot  of  the  A-shock  structure  is  accurately  predicted  by  the 
computation,  the  large  gradients  in  the  compression  waves  are  smoothed  out.  The  compu¬ 
tation  shows  that  both  the  front  and  the  rear  foot  of  the  A-shock  are  approximately  planar 
within  the  accuracy  of  the  numerical  interpolation.  The  two  systems  of  the  compression 
waves  merge  at  about  y  —  2.0fao. 

The  measured  and  computed  surface  streamlines  are  shown  in  Figs.  13  and  14  respec¬ 
tively.  As  in  the  iin  case,  all  of  the  qualitative  features  of  the  surface  flow  visualization  are 
predicted  by  the  computation.  The  computed  upstream  influence  and  line  of  coalescence  are 
in  reasonable  agreement  with  experiment  with  relative  underprediction  of  about  10%.  The 
divergence  line  (attachment)  is  located  on  the  model  surface.  The  analysis  of  the  surface 
pattern  is  similar  to  that  provided  for  the  sharp  fln  above. 

Computed  mean  streamlines  are  displayed  in  Figs.  15  and  16.  In  Fig.  15,  three  sets  of 
streamlines  are  shown,  the  first  released  within  the  computational  sublayer  near  the  wall,  the 
second  and  third  sets  released  at  y  =  0.5^ao  V  —  ^-26oo  respectively.  The  overall  flow 
structure  is  somewhat  similar  to  that  for  the  sharp  fln.  The  first  set  of  particles  form  the 
line  of  coalescence,  the  second  set  is  entrained  into  the  vortical  structiue  while  the  the  third 
set  is  lifted  by  the  model  and  moves  downstream  along  the  comer  surface.  Fig.  16  displays 
a  side  view  of  the  previous  plot  with  only  the  second  and  third  sets.  Other  sets  of  trues  are 
consistent  with  the  above  results. 

11*2. S  Comments  on  Quasiconical  Similarity 

For  the  purposes  of  discussion,  the  surface  pressure  of  the  two  models  are  compared  in  Fig.  17. 
The  experimental  data  is  superimposed.  The  experiment  clearly  shows  that  the  upstream 
portion  of  the  surface  pressure  is  almost  identical  for  the  interaction  generated  by  the  two 
models.  The  differences  caused  by  the  particular  geometry  of  the  model  appear  only  behind 
the  inviscid  shock  wave.  The  computed  STirface  pressure  of  the  two  interactions  also  shows 
some  similarity  upstream  of  the  shock  wave. 

The  flowfield  structure  is  investigated  using  the  particle  tracing  technique  and  MO  VIE.BYU 
color  graphics  software.  In  Fig.  18,  two  flow  ribbons  and  two  flow  stream  surfaces  of  the  17.5 
deg  fin  interaction  flowfield  are  displayed.  The  red  surface  is  the  ’surface  of  separation’  which 
emanates  from  the  line  of  coalescence  and  rises  from  the  flat  plate  into  the  vortical  core.  The 
blue  surface  is  called  ’surface  of  attachment’  which  intersects  the  flat  plate  at  the  line  of  di¬ 
vergence.  The  lower  ribbon  begins  to  turn  near  the  model  apex  and  moves  below  the  surface 
of  separation  twisting  along  the  way.  The  upper  flow  ribbon  which  is  originally  above  the  sur¬ 
face  of  attachment  is  not  entrsdned  into  the  vortical  structure  and  moves  downstream  almost 
parallel  to  the  fin  siuface.  The  overall  flowfield  structme  is  apparent.  Particles  originating 
between  the  plate  and  the  surface  of  attachment  are  entrained  into  the  vortical  structure 
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while  those  originating  above  the  surface  of  attachment  flow  over  the  vortical  structure. 

II-3  Research  Task  No.  3:  Investigate  Methods  of  Control  and  Modifi¬ 
cation  of  3-D  Turbulent  Interactions 

As  indicated  previously,  the  motivation  for  this  research  task  lies  in  the  existence  of  an 
undesirable  large  vortical  structure  in  3-D  turbulent  interactions.  Several  techniques  may  be 
appropriate  for  flow  control  including  1)  bleed  (suction)  to  remove  low  energy  fluid  in  the 
boundary  layer,  2)  generation  of  vorticity  of  opposite  sign  (e.g.,  with  the  symmetric  double 
wedge  configuration  (Mee  and  Stalker  [18])),  3)  blowing  (normal  or  wall- jet)  to  enhance  the 
kinetic  energy  of  particles  in  the  boundary  layer,  and  4)  geometrical  modifications  such  as 
raising  the  sharp  fin  a  finite  distance  above  the  surface^. 

The  focus  the  present  research  is  two-fold: 

•  Bleed  (suction):  Bleed  has  been  traditionally  employed  in  high  speed  aircraft  inlets 
to  prevent  boundary  layer  separation.  Conventional  theory  dictates  that  the  effect  of 
suction  is  to  remove  the  low  speed  fluid  in  the  boimdary  layer  before  it  separates  from 
the  surface. 

•  Shock-shock  Intersections:  Mee  and  Striker  [18]  have  concluded  from  experimental 
observations  that  intersecting  shocks  can  produce  a  given  overall  pressure  rise  with  lesser 
likelihood  of  separation  than  an  equivalent  strength  single  shock  interaction.  Their 
experiments  were  restricted  to  weak  interactions.  This  concept  is  unproven,  however, 
and  additional  computational  and  experimental  research  is  required  to  determine  its 
validity. 

II-S.l  EflTeet  of  Bleed 

A  series  of  computations  on  the  20**  single  fin  configuration  (SFC)  with  Mach  and  Re3molds 
numbers  3  and  9  x  10^,  with  porous  suction  over  a  triangular  region  extending  from  the 
line  of  upstream  infiuence  to  the  theoretical  inviscid  shock  location  (Region  I  in  Fig.  19)  has 
been  described  in  previous  literature  [19].  The  maTiTniiTn  bleed  corresponded  to  removal  of 
approximately  25  %  of  the  upstream  boundary  layer.  The  computed  streamline  structure 
was  observed  to  be  essentially  independent  of  the  stirface  bleed  flow  rate,  although  the  line 
of  coalescence  moved  downstream  toward  the  theoretical  inviscid  shock  footprint  (TISF) 
and  upstream  influence  was  reduced.  The  effect  of  bleed  was  focused  principally  in  a  small 
fraction  of  the  boundary  layer.  This  represents  a  significant  and  surprising  restilt  and  brings 
into  question  the  efficacy  of  boundary  layer  bleed  for  sidewall  interactions  in  high  speed 
inlets.  Numerical  simulation  of  bleed  on  a  different  area  -  the  region  on  the  plate  between 
the  TISF  and  the  base  of  the  fin  (Region  II  in  Fig.  19)  and  performed  in  FY88  are  described. 
Approximately  8%  of  the  incoming  boundary  layer  is  bled  off  from  the  region  on  the  flat  plate 
extending  between  the  theoretical  inviscid  shock  line  and  the  surface  of  the  sharp  fin.  The 
influence  of  bleed  in  this  region  on  the  vortical  structure  is  also  found  to  be  insignificant.  A 
detailed  description  of  the  results  may  be  found  in  Gaitonde  and  Knight  [20]*. 

^S.  BogdonofT,  Private  communication 
‘See  Section  VII 
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The  recent  experiment  of  Barnhart  et  al  [21]  with  suction  on  a  roughly  rectangular  area 
covering  a  strip  on  either  side  of  the  theoretical  shock  footprint  (Fig.  20)  on  an  incipiently 
separated  configuration  (8°  SFC,  Mach  3,  RrSo,,  =  ^  10^)  is  extensive  enough  to 

provide  general  guidelines.  Nevertheless,  their  observation  of  beneficial  effects  of  suction 
warrants  numerical  investigation.  This  computation  is  complete  and  results  are  described 
in  [20].  Briefly,  comparison  with  experimental  data  and  analysis  of  computed  results  indicates 
that  the  effect  of  bleed  is  primarily  to  reduce  the  surface  angularity  and  upstream  influence. 
The  major  portion  of  the  flow  is  not  influenced  by  suction. 

The  overall  effect  of  suction,  based  on  the  above  computations  on  the  strong  (20  deg)  and 
weak  (8  deg)  configurations,  may  be  summarized  as  follows: 

•  Effect  on  s\irface  streamlines:  For  the  stronger  interaction,  the  effect  of  bleed  upstream 
of  the  shock  footprint  (Region  I)  is  to  reduce  the  surface  angularity  of  the  flow.  The 
line  of  coalescence  observed  in  the  absence  of  bleed  [7]  persists  although  its  angle  with 
the  shock  footprint  is  significantly  reduced.  This  reduction  is  proportional  to  bleed 
magnitude.  Region  il  bleed  has  no  effect  on  surface  streamlines  which  resemble  those 
in  the  absence  of  suction.  For  the  weaker  interaction,  surface  yaw  angles  also  display 
lower  values  in  presence  of  suction.  A  general  statement  about  the  effect  of  bleed 
on  separation  cannot  be  made  for  the  weaker  configuration  since  it  is  only  incipiently 
separated  even  in  the  absence  of  bleed. 

•  Effect  on  surface  variables:  Region  I  bleed  results  in  lower  upstream  influence  and 
higher  rate  of  pressure  rise  in  the  interaction  although  the  total  pressure  rise  remains 
the  same.  Region  II  bleed  does  not  display  any  effect  on  surface  pressure  rise.  For 
the  weak  interaction,  the  results  are  similar  to  those  for  Region  I  bleed  on  the  strong 
interaction.  In  the  absence  of  further  computation,  it  is  not  possible  to  distingtiish 
between  the  effects  of  upstream  and  downstream  (of  the  shock  footprint)  bleed.  The 
effect  of  all  regions  of  bleed  considered  (on  both  the  strong  and  the  weak  configuration) 
is  to  increase  local  skin  friction  drastically.  It  is  evident  that  a  significant  drag  penalty 
is  associated  with  bleed. 

•  Effect  on  shock  structure:  Region  I  bleed  results  in  a  tightening  of  the  upstream  leg 
of  the  A-shock  structure  as  a  result  of  reduced  displacement  thickness  of  the  interac¬ 
tion.  This  reduction  is  evident  also  by  the  elimination  of  the  pitot  pressure  overshoot 
observed  in  the  interaction  in  the  absence  of  bleed.  Region  II  bleed  has  no  effect  on 
either  the  upstream  or  the  downstream  compression  fans.  For  the  weak  interaction,  the 
upstream  compression  fan  is  seen  to  disappear  completely.  This  is  probably  due  to  a 
combination  of  the  relatively  high  values  of  bleed  utilized  and  the  inherent  weakness  of 
the  interaction. 

•  Effect  on  vortical  structure:  Neither  bleed  region  has  any  significant  effect  on  the 
vortical  structure  observed  in  the  strong  interaction  although  larger  portions  of  the 
incoming  boundary  layer  are  ingested  with  increased  bleed  magnitude.  Only  particles 
whose  no  bleed  trajectories  pass  close  to  the  bleed  region  are  ingested.  Although  no 
statement  can  be  made  about  the  effect  of  bleed  on  the  vortical  structure  for  the  weaker 
interaction,  it  is  observed,  that  particles  escaping  the  influence  of  bleed  follow  paths 
similar  to  those  in  the  no  bleed  configuration. 
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II-S.2  Double  shock  interactions 

Following  the  conclusions  of  Mee  and  Stalker  that  shock-shock  intersections  may  be  capa¬ 
ble  of  producing  a  given  pressure  rise  with  less  likelihood  of  separation  than  a  single  shock 
interaction,  a  systematic  study  of  symmetric  double  fin  configurations  is  initiated.  A  calcula¬ 
tion  was  performed  previously  [19]  for  the  strongest  intersecting  shock  configuration  of  Mee 
and  Stalker  corresponding  to  incoming  Mach  and  Reynolds  numbers  1.85  and  7.8  x  10^  and 
symmetric  sharp  fins  with  angles  of  5  deg  each.  The  results  indicated  overall  good  agreement 
with  experiment  with  one  interesting  anomaly  that  the  computation  overpredicted  the  extent 
of  upstream  influence  by  approximately  13%.  This  is  in  contrast  to  previous  3-D  turbulent 
interaction  studies  [7]  where  computations  have  consistently  underpredicted  upstream  influ¬ 
ence.  A  more  systematic  study  of  stronger  interactions  is  described  in  this  research.  The 
Mach  and  Reynolds  numbers  are  2.95  and  2  x  10^  respectively  for  S3rmmetric  fin  angles  of 
4  deg  and  8  deg.  The  flow  parameters  are  chosen  after  consideration  of  experimental  facilities 
available  at  the  Gas  Dynamics  Laboratory  at  Priceton  University.  The  4x4  configuration 
is  a  completely  weak  interaction  for  the  domain  under  consideration.  In  contrast,  the  8x8 
configuration  displays  a  surface  flow  pattern  with  yaw  angles  larger  than  the  primary  shock 
angle  and  an  interesting  node-like  structure  (whose  validity  is  suspect)  near  the  trailing  edge 
of  the  domain.  The  consequent  separated  flow  is  examined  with  particle  traces  (Section  VII). 
Particles  released  close  to  the  plate  but  outside  the  sublayer  follow  the  general  pattern  of  the 
surface  streamlines  up  to  the  node  at  which  point  they  rise  and  move  parallel  to  the  center 
line.  Particles  released  in  the  vicinity  of  the  line  of  symmetry  at  various  heights  move  parallel 
to  it  (as  expected  from  the  imposed  boundary  conditions)  but  rise  up  and  away  from  the  plate 
to  form  an  arch-like  trace  at  the  downstream  end  of  the  domain.  Particles  released  near  the 
comer  formed  by  the  fin  and  the  plate  pitch  toward  the  plate  and  yaw  toward  the  line  of 
symmetry  straightening  out  as  they  approach  it.  At  higher  distances  from  the  plate,  particles 
released  near  the  fin  pitch  toward  the  plate  over  the  particles  beneath  yawing  toward  the  cen¬ 
terline.  The  arch-like  structure  formed  by  the  particles  released  near  the  centerline  broadens 
so  that  the  overall  flow  approaches  the  boundary  layer  flow  (due  to  the  fin)  expected  at  large 
distances  from  the  plate.  It  is  anticipated  that  equivalent  experiments  will  be  performed  in 
the  near  future. 

II-4  Research  Task  No.  4:  Examination  of  inviscid  and  viscous  effects 
in  3-D  turbulent  interactions 

The  flowfleld  results  luing  the  two  different  turbulence  models  for  both  fin  and  swept  comer 
configurations  have  been  found  to  be  in  close  agreement  except  for  modest  differences  in  a 
small  fraction  of  the  boundary  layer  adjacent  to  the  siuface.  This  relative  insensitivity  of 
the  flowfleld  to  the  type  of  turbulence  model  employed  in  the  computations  suggests  that  the 
interaction  is  primarily  an  inviscid  and  rotational  phenomenon,  with  viscous  influence  being 
limited  to  regions  very  close  to  the  wall  or  surface.  The  present  study  was  initiated  to  verify 
this  hypothesis  throug;-  extensive  investigation  of  the  sharp  fin  and  swept  compression 
comer  flowflelds  computed  using  the  two  different  turbulence  models. 

The  analysis  of  the  flow  past  the  20  deg  fln  and  the  (24,40)  swept  compression  comer 
both  at  Mach  3,  Re^  ~  9  x  10®  is  complete.  The  actual  flows  were  computed  separately  by 
Knight  and  Horstman  who  utilize  the  Baldwin- Lomax  and  Jones-Launder  turbulence  models 


12 


respectively. 

The  following  approach  is  adopted  in  the  analysis  of  the  flowhelds: 

1.  A  niunerical  code  is  developed  to  integrate  the  mean  kinetic  energy  equation  along 
streamlines  in  the  flow  and  to  compute  the  contributioiu  to  the  rate  of  change  of 
mean  kinetic  energy  by  viscous  and  inviscid  terms.  Using  this  code,  several  streamlines 
originating  at  various  y^locations  within  and  outside  the  boundary  layer  are  examined. 
The  results  for  both  sharp  fln  and  swept  compression  comer  flows  indicate  that  although 
inviscid  effects  in  general  govern  the  dynaxnics  of  the  interaction  along  the  streamlines, 
viscous  influence  remains  significant  even  at  points  considerably  distant  from  the  plate 
or  surface. 

2.  A  ‘planewise’  investigation  of  the  interaction  flowfields  is  carried  out  in  order  to  directly 
determine  the  relative  importance  of  viscous  and  inviscid  contributions  in  various  re¬ 
gions  of  the  flow.  The  magnitudes  of  the  viscous  and  inviscid  terms  (of  the  mean  kinetic 
energy  equation)  are  compared  directly  at  points  lying  on  various  y  —  z  planes  (for  the 
sharp  fin)  and  at  points  lying  on  various  z  —  y  planes  (for  the  swept  compression  comer). 
It  is  observed  that  although  the  interaction  is  governed  primarily  by  inviscid  effects, 
viscous  influence  is  significant  even  in  regions  of  the  flow  distant  from  the  surface. 

Color  graphics  is  extensively  utilized  to  visualize  the  results  of  this  ‘planewise’  analysis. 
Figs.  21  and  22  display  the  magnitudes  of  the  inviscid  and  viscous  terms  respectively  for 
the  sharp  fin  interaction  (a  =  20  deg,  Mach  3,  Res  =  9  x  10^),  computed  with  the  the 
Baldwin-Lomax  turbulence  model,  on  a  y  —  z  plane  (normal  to  the  streamwise  direction) 
and  at  a  location  z  —  IS^eo'  The  inviscid  (pressure)  influence  is  concentrated  towards  the 
fleestream  side  of  the  interaction  (i.e.  away  from  the  fin  surface).  The  viscous  effects  are 
concentrated  to  the  immediate  left  of  the  region  of  inviscid  influence.  Similar  results  have 
been  observed  for  other  y  —  z  planes  at  various  z  locations.  The  flowfield  computed  using 
the  Jones-Launder  model  shows  similar  results  with  regard  to  the  concentration  of  inviscid 
influence  in  the  flow  but  viscous  influence  is  observed  to  be  concentrated  in  a  fraction  of  the 
boundary  layer  adjacent  to  the  surfaces  of  the  fin  and  the  plate. 

The  results  for  the  swept  comer  interaction  computed  using  the  Baldwin-Lomax  model 
parallel  those  for  the  sharp  fin  for  both  viscous  and  inviscid  terms.  The  flowfield  computed 
using  the  Jones-laimder  model  for  this  interaction  agrees  as  far  as  the  inviscid  terms  are 
concerned  but  differs  with  respect  to  the  region  of  concentration  of  viscom  influence. 

A  more  complete  discussion  of  the  approach  to  the  problem  and  the  results  of  this  research 
task  is  included  in  the  attached  paper  titled  “Viscous/Inviscid  Influence  in  3D  Shock  Wave- 
Turbulent  Boundary  Layer  Interactions”,  submitted  for  the  AlAA  Meeting  in  January,  1989 
and  included  in  Section  VII. 
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Ill  Research  Program  for  FY89  (1  October  1988  -  30  Septem> 
ber  1989) 

in-l  Research  Task  No.  1:  Develop  and  Validate  Theoretical  Models 
for  3-D  Shock  Wave-Turbulent  Boundary  Layer  Interactions 

A  series  of  computations  will  be  performed  during  FY89  for  the  3>D  sharp  iin  configuration 
at  Mach  7  at  fixed  Reynolds  number  for  several  fin  an^es  a.  Additional  computations  may 
be  performed  at  Mach  10.  The  flow  conditions  will  be  selected  to  agree  with  the  planned 
experiments  at  NASA  Ames.  Although  the  precise  requirements  for  the  number  of  grid  points 
and  CPU  time  depends  on  the  specific  case,  experience  has  indicated  that  typically  100,000 
grid  points  and  10  CYBER  205  CPU  hours  are  required  [7]  for  each  case. 

The  computed  results  will  be  compared  with  experimental  data.  Present  plans  include 
measurements  of  surface  pressure  and  heat  transfer,  surface  oil  film  visualisaticm,  and  bound¬ 
ary  layer  profiles  of  pitot  and  static  pressure,  velocity  and  turbulence  fluctuations.  The  data 
and  graphics  postprocessing  will  be  performed  at  the  Rutgers  College  of  Engineering  Super¬ 
computer  Remote  Access  and  Graphics  Facility. 

in-2  Research  Task  No.  2:  Investigate  Methods  for  Control  and  Modi¬ 
fication  of  3-D  Turbulent  Interactions 

The  proposed  research  activity  is  strongly  dependent  on  the  avulability  of  experimental  data 
and  win  focus  on  the  double  fin  configuration  and  the  lifted  fin  configuration  (consisting  of  a 
fin  slightly  separated  from  the  flat  plate  thus  permitting  flow  bleed  through  the  comer).  The 
computations  on  the  double-fin  configuration  utilized  parameters  chosen  in  anticipation  of 
experiments  to  be  performed  at  the  Princeton  Gas  Dynamics  Laboratory.  In  particular,  the 
lengths  of  the  fins  and  the  spanwise  separation  (Fig.  3)  were  tentatively  based  on  dimensions 
of  the  wind  tunnel  to  be  employed.  The  experimental  apparatus  to  be  actually  utilized  is 
currently  in  advanced  stages  of  calibration  [22].  The  apparatus  differs  from  the  configuration 
employed  in  two  ways; 

•  Model  Structure:  The  experimental  setup  exhibits  both  vertical  and  horizontal  planes  of 
symmetry  thus  resembling  more  closely  an  inlet  configuration.  Horizontal  symmetry  is 
achieved  by  utilizing  a  top  plate  symmetric  with  the  bottom  plate.  The  physical  domain 
employed  for  double  shock  interactions  in  Section  II-3  does  not  possess  a  horizontal 
plane  of  symmetry  and,  consequently  asymptotic  2-D  boundary  conditions  are  enforced 
at  the  top  boimdary. 

•  Dimensions:  The  actual  dimensions  of  the  model  are  currently  being  redetermined 
based  on  practical  wind-tuimel  operation  constraunts. 

The  systematic  series  of  experiments  initiated  in  FY87  will  be  continued  in  FY89  with  close 
collaboration  with  Prof.  BogdonofF  of  Princeton  University. 

A  preliminary  experimental  investigation  into  the  effect  of  comer  bleed  [23]  has  been 
recently  completed  at  the  Princeton  Gas  Dynamics  Laboratory.^.  Although  lifting  the  fin  off 

*S.  Bogdonoff,  Private  coininunication. 
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the  plate  does  not  alter  the  components  of  the  physical  model  (iin  and  plate)  in  any  drastic 
fashion,  significant  complexities  in  grid  generation  and  boundary  conditions  (e.g.,  under  the 
fin)  may  arise.  An  initial  investigation  into  the  feasibility  of  such  a  computation  will  be 
carried  out. 

in-3  Research  Task  No.  3:  Examination  of  Inviscid  and  Viscous  Effects 
in  3-D  Turbulent  Interactions 

Efforts  will  continue  on  the  quantitative  analysis  of  inviscid  and  viscous  effects  for  the  3-D 
swept  compression  comer  and  the  3-D  sharp  fin  flowfields.  These  activities  will  include: 

•  Examination  of  the  error  in  the  evaluation  of  the  mean  kinetic  energy  for  the  swept 
compression  comer  for  (a,  A)  =  (24,40)  deg  at  Reynolds  number  2.6  x  10^. 

•  The  flowfields  computed  to  investigate  the  quasiconical  free  interaction  principle  (Sec¬ 
tion  n-2)  will  be  analysed  for  viscous-inviscid  effects 

•  Other  flow  equations  (e.g.,  the  mean  vorticity  equation)  will  be  examined 

High  resolution  color  graphics  visualization  will  be  employed  for  presentation  of  the  re¬ 
sults. 
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Figure  2:  Flowfield  for  swept  comer 
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Figure  3:  Double  fin  configuration 


Figure  4:  Location  of  experimental  surveys  —  17.5  deg  fin 
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Figure  6  (a)  (1  of  3):  Pitot  pressure  at  Station  12  -  17.5  deg  fin 
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Figure  6:  Static  pressure  comparison  at  Station  12  -  17.5  deg  lin  (3  of  3) 
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Figure  11:  Surface  pressure  at  z,  =  lOfgo  for  (30,60)  comer 
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Figure  17:  Surface  pressure  profiles  for  17.5<fep  Fin  and  (30,60)  Swept  Compression  Comer 
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Theoretical  Inviscid  Shock 


Figure  19;  Bleed  Regions  —  20  deg  Single  Fin  Configuration  (SFC) 
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Figure  20:  Bleed  Regions  —  8  deg  Single  Fin  Configuration  (SFC) 
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Numerical  Investigation  of  Some  Control  Methods  for  3-D  Turbulent 

Interactions  due  to  Sharp  Fins 


Gaitonde^ 

Doyle  Knight’ 

Department  of  Mechanical  and  Aerospace  Engineering 
Rutgers  University 
P.O.  909,  MAE 
Piscataway,  NJ  08855 


Abstract 

This  paper  describes  continuing  numerical  re¬ 
search  efforts  on  the  modification  and  control  of  3-D 
shock  wave-turbulent  boundary  layer  interactions  due 
to  supersonic  flow  past  a  sharp  fin  mounted  on  a  flat 
plate.  In  supplement  to  previous  numerical  experiments 
(20*  fin,  Mach  3,  Rtf  =  9  x  10‘),  the  effect  of  bleeding 
off  roughly  5%  k  8%  of  the  incoming  boundary  layer 
from  two  distinct  bleed  regions  is  studied.  Region  I 
covers  the  area  on  the  flat  plate  between  the  line  of 
upstream  influence  (as  determined  by  computed  initial 
pressure  rise)  and  the  theoretical  inviscid  shock  foot¬ 
print  (TISF)  while  Region  II  covers  the  remaining  por¬ 
tion  on  the  plate  between  the  TISF  and  ti  e  base  of  the 
fin.  Numerical  pitot  pressure  and  yaw  angles  are  com¬ 
pared  with  experimental  data  obtained  in  the  absence 
of  bleed.  The  results  indicate  that  the  influence  of  both 
suction  schedules  on  the  overall  flowfield  and  the  previ¬ 
ously  observed  vortical  structure  is  remarkably  modest 
with  significant  effects  restricted  to  the  region  a((jacent 
to  the  plate  near  the  bleed  region.  In  a  separate  com¬ 
putation,  for  an  incipiently  separated  configuration  (8* 
fin,  Mach  3,  Rcf  =  8.5  x  lO'),  23%  suction  is  applied 
on  an  area  across  the  shock  corresponding  to  the  ex¬ 
perimental  configuration  of  Barnhart  et  aL  Comparison 
with  experimental  data  and  analysis  of  computed  results 
indicates  that  the  effect  of  bleed  is  primarily  to  reduce 
surface  angularity  and  upstream  influence.  The  major 
portion  of  the  flow  is  not  influenced  by  suction. 

Following  the  conclusions  of  Mee  and  Stalker  that 
shock-shock  intersections  may  be  capable  of  producing 
a  given  pressure  rise  with  less  likelihood  of  separation, 
a  systematic  study  of  symmetric  double  fin  configura¬ 
tions  is  initiated.  Two  configurations  with  fin  angles 
4  deg  X  4  deg  and  8  deg  x  8  deg  respectively  at  Mach  2.95 
and  Reynolds  number  2  x  10’  are  simulated.  Results 
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are  presented  in  the  form  of  surface  streamlines,  surface 
pressures  and  particle  traces.  The  4x4  configuration 
is  a  completely  weak  interaction  for  the  domain  under 
consideration.  In  contrast,  the  8x8  configuration  dis¬ 
plays  surface  flow  turning  angles  larger  than  the  primary 
shock  angle.  The  consequent  separated  flow  is  examined 
with  particle  traces. 

1  Introduction 

Research  on  3-D  shock  wave  -  turbulent  boundary 
layer  interactions  (‘‘3-D  turbulent  interactions”)  finds 
utility  in  a  number  of  high  speed  eterodynamic  applica¬ 
tions  such  as  aircraft  components.  A  supersonic  inlet 
flowfield,  for  example,  is  characterised  by  a  pattern  of 
oblique  shocks  formed  by  the  general  curvilinear  shape 
of  the  tamp  and  cowl  surfaces,  that  interact  with  turbu¬ 
lent  boundary  layers  on  the  walls.  Dimensionless  geome¬ 
tries  such  as  the  3-D  sharp  fin  (as  in  this  research)  and 
the  3-D  swept  compression  corner  are  often  employed  to 
isolate  the  physics  of  such  interactions  from  the  geomet¬ 
rical  complexity.  The  principal  parameters  for  the  sharp 
fin  (Fig.  1)  configuration  are  the  Mach  number  Mgo  ,  the 
Reynolds  number,  Rcf^  based  upon  the  boundary  layer 
thickness  Soo  the  leading  edge  of  the  fin,  the  thermal 
boundary  conditions  and  the  fin  angle  a. 

Progress  has  been  achieved  in  analysing  3-D  tur¬ 
bulent  interactions  with  analytic,  experimental  and  nu¬ 
merical  techniques  applied  separately  or  in  coqjunction 
with  each  other.  Analytical  investigations  have  been 
restricted  primarily  to  weak  interactions  in  which,  by 
definition,  there  is  no  separation  [1,  2,  3,  4].  Such  stud¬ 
ies  have  provided  valuable  insight  into  the  various  scales 
involved  in  complex  interactions  resulting  in  better  un¬ 
derstanding  of  required  resolution  for  numerical  and  ex¬ 
perimental  effort. 

Experimental  results  in  the  form  of  surface  pres¬ 
sures  and  surface  flow  visualisations  have  been  available 
for  over  two  decades[5,  6,  7,  8,  9,  10].  In  recent  years, 
improvements  in  instrumentation  have  made  it  possible 
to  meuure  flow  variables  (such  as  pitot  pressure  and 
yaw  angles)  in  the  boundary  layer  as  well  [11,  12,  13, 
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14].  Meuurements  of  skin  friction  by  laser  interferom¬ 
etry  [15]  have  helped  compare  the  accuracy  of  different 
eddy  viscosity  models.  In  addition  to  providing  results 
of  practical  value,  experimental  observations  provide  a 
valuable  benchmark  for  code  validation. 

The  objective  of  current  research  efforts  is  to 
achieve  the  ability  to  not  only  predict  variables  of  prac¬ 
tical  interest  (surface  pressures,  sldn  friction  and  heat 
transfer  rates),  but  also  overall  flow  characteristics  such 
as  upstream  influence,  symmetricity  (conical  or  cylin¬ 
drical  [16]),  existence  of  separation  and  flow  structures 
on  various  scales.  With  the  advent  of  relatively  eas¬ 
ily  accessible  high-speed  computers,  computational  re¬ 
search  efforts  have  become  a  very  cost-efficient  means  to 
"reduce  the  required  [experimental]  test  matrix  to  the 
smallest  number  of  configurations  [17]”. 

Computations  [18,  19,  20,  21,  22]  typically  employ 
the  Reynolds  averaged  Navier-Stokes  equations  coupled 
with  an  algebraic  [23]  or  two-equation  [24]  turbulence 
model.  Numerical  results  for  the  sharp  fin  [14]  and  the 
swept  compression  corner  [25,  26]  have  shown  overall 
good  agreement  with  experimental  measurements,  with 
modest  local  discrepancies  in  the  prediction  of  upstream 
influence  and  skin  fnction  values.  Analysis  of  the  com¬ 
puted  flowfields  confirm  the  principal  flow  feature  —  a 
large  vortical  structure  aligned  with  the  cornet  in  agree¬ 
ment  with  the  flowfield  models  of  Token  [27]  and  Kubota 
and  Stollery  [8].  A  three  dimensional  surface  of  separao 
tion  (Fig.  2)  emanates  from  the  line  of  coalescence  (sep¬ 
aration),  and  spirals  into  the  vortical  center.  A  second 
surface,  emanating  from  upstream,  intersects  the  wall 
at  the  line  of  divergence  (attachment),  and  defines  the 
extent  of  the  fluid  entrained  into  the  vortical  structure. 

Such  a  flowfield  -  specifically  the  large  vortical 
structure  and  separated  flow  -  is  undesirable  in  applica¬ 
tions  such  as  aircraft  inlets  whose  function  is  to  provide 
a  nearly  uniform  subsonic  flow  with  high  total  pressure 
recovery  at  the  compressor  face.  The  objective  of  this 
research  work  is  to  study  possible  perturbations  to  the 
boundary  conditions  which  might  influence  the  flowfield 
in  a  positive  fashion  and  thus  provide  control.  The  need 
for  boundary  layer  control  studies  has  been  recognized 
in  the  literature; 

Despite  the  obvious  importance  of  boundary- 
layer  control  in  high-speed  air-breathing 
propulsion  inlets,  only  two  basic  experimental 
studies  of  this  problem  involving  swept  inter¬ 
actions  have  been  found  in  the  literature  both 
employing  sharp  fin  shock  generators  [28]. 

Several  possible  methods  of  control  may  be  pro¬ 
posed: 

•  Bleed  (suction):  Boundary  layer  suction  serves  to 
remove  low  speed  fluid  in  the  boundary  layer  before 


it  separates  from  the  surface.  In  inlets,  suction  is 
traditionally  employed  on  the  tamp,  cowl  and  side- 
walls  to  prevent  flow  separation  in  the  vicinity  of 
the  intersection  of  the  shock  waves  and  walls. 

•  Blowing;  Experimental  results  [12]  at  several  Mach 
numbers  mdicate  that  tangential  blowing  is  benefi¬ 
cial  in  retarding  separation.  The  effect  of  blowing 
is  to  enhance  the  energy  of  the  retarded  boundary 
layer  fluid. 

•  Shock-shock  Interaction:  Mee  and  Stalker  [4]  have 
theorised  from  an  experimental  study  of  weak 
shock-shock  interactions  with  a  turbulent  bound¬ 
ary  layer  that  intersecting  shocks  can  produce  a 
given  overall  pressure  rise  with  less  likelihood  of 
separation  than  an  equivalent  strength  single  shock 
interaction. 

•  Modification  of  Geometrsr;  It  has  been  suggested 
that  minor  modifications  in  the  3-D  configuration 
itself  may  exert  significant  influence.  One  pro¬ 
posed  [29]  envisages  lifting  the  fin  off  the  surface 
of  the  flat  plate  effectively  permitting  the  flow  to 
bleed  off  &om  the  comer  formed  by  the  fin  and  the 
plate. 

•  Vorticity  Modifications  in  Incoming  Flow:  Since  the 
principal  flow  feature  is  a  vortical  structure,  it  may 
be  possible  to  apply  control  through  the  introduc¬ 
tion  of  additional  vorticity  in  the  longitudenal  or 
spanwise  direction. 

This  paper  describes  a  continuing  research  ef¬ 
fort  focused  on  a  numerical  investigation  of  the  first 
and  third  of  these  options.  For  suction,  the  single  fia 
configuration  (SFC  -  Fig.  1)  is  utilised.  A  series  of 
computations  on  the  20"  SFC,  with  porous  suction  over 
a  triangular  region  extending  &om  the  line  of  upstream 
influence  to  the  theoretical  inviscid  shock  location  (Re¬ 
gion  I  in  Fig.  3)  has  been  described  in  previous  litera¬ 
ture  [30].  The  maximum  bleed  corresponded  to  removal 
of  approximately  25  %  of  the  upstream  boundary  layer. 
The  computed  streamline  structure  was  observed  to  be 
essentially  independent  of  the  surface  bleed  flow  rate, 
although  the  line  of  coalescence  moved  downstream  to¬ 
ward  the  TISF  and  upstream  influence  was  reduced. 
The  effect  of  bleed  was  focused  principally  in  a  small 
fraction  of  the  boundary  layer.  This  represents  a  signif¬ 
icant  and  surprising  result  and  brings  into  question  the 
efficacy  of  boundary  layer  bleed  for  sidewall  interactions 
in  high  speed  inlets.  Numerical  results  of  the  applica¬ 
tion  of  bleed  on  a  different  area  -  the  region  on  the  plate 
between  the  TISF  and  the  base  of  the  fin  (Region  II  in 
Fig.  3)  are  described  in  this  paper  (Section  6.1).  The 
recent  experiment  of  Barnhart  et  al  [31]  with  suction  on 
an  incipiently  separated  configuration  (8"  SFC,  Mach  3, 
Res^  =  6.5  X  10®)  is  not  extensive  enough  to  provide 
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geneial  guidelenes.  Nevertheless,  their  observation  of 
beneficial  effects  of  suction  warrants  numerical  investi¬ 
gation  (Section  6.2). 

The  second  set  of  computations  deals  with  crossed 
shock  wave  -  turbulent  boundary  layer  interactions. 
The  concept  of  employing  intersecting  shocks  (produced 
by  the  symmetric  double  jin  configuration  (SDFC)  in 
Fig.  4)  is  based  on  the  experimental  observation  that 
the  limiting  streamline  deflection  angles  downstream  of 
the  reflected  shocks  are  less  tham  the  angles  downstream 
of  the  primary  shock  [32].  This  concept  is  unproven, 
however,  and  additional  computational  and  experimen¬ 
tal  research  is  required  to  determine  its  validity.  It  is 
emphasized  that  suction  in  not  imposed  in  the  double 
fin  configurations.  A  calculation  was  performed  previ¬ 
ously  [30]  for  the  strongest  intersecting  shock  configura¬ 
tion  of  Mee  and  Stalker  corresponding  to  incoming  Mach 
and  Reynolds  numbers  1.85  and  7.8  x  10^  and  symmetric 
sharp  fins  with  angles  of  5  deg  each.  The  results  indi¬ 
cated  overall  good  agreement  with  experiment  with  one 
interesting  anomaly  that  the  computation  overpredicted 
the  extent  of  upstream  influence  by  approximately  13%. 
This  is  in  contrast  to  previous  3-D  turbulent  interac¬ 
tion  studies  [14]  where  computations  have  consistently 
underpredicted  upstream  influence.  A  more  systematic 
study  of  stronger  interactions  is  described  in  this  paper. 
The  Mach  and  Reynolds  numbers  are  2.9S  and  2  x  10* 
respectively  for  symmetric  fin  angles  of  4  deg  and  8  deg. 
It  is  anticipated  that  equivalent  experiments  will  be  per¬ 
formed  at  the  Princeton  Gas  Dynamics  Laboratory  in 
the  future. 

2  Governing  Equations  and  Numerical  Model 

The  governing  equations  are  the  full  mean  com¬ 
pressible  Navier-Stokes  equations  using  mass  averaged 
variables  and  strong  conservation  form  [33].  The  molec¬ 
ular  and  turbulent  Prandtl  numbers  are  0.73  and  0.9  re¬ 
spectively.  Turbulence  is  modeled  through  the  inclusion 
of  the  two-layer  algebraic  turbulent  eddy  viscosity  model 
of  Baldwin  and  Lomax  [23]  with  the  mixing  length  (/) 
as  specified  by  the  formula  of  Buleev  [34].  The  effect  of 
bleed  is  incorporated  through  the  bleed  correction  fac¬ 
tor  of  Cebeci  [35]  u  described  in  [30].  This  turbulence 
model  has  been  validated  for  a  variety  of  boundary  layer 
flows  [36]  with  zero,  favorable  and  adverse  pressure  gra¬ 
dients  in  the  presence  of  suction. 

A  3-D  coordinate  transformation  is  employed  as 
in  [14]  to  facilitate  the  application  of  the  numerical 
algorithm  and  the  boundary  conditions.  The  hybrid 
explicit-implicit  algorithm  of  Knight  [33]  is  employed 
to  solve  the  governing  equations.  The  implicit  algo¬ 
rithm  employs  Keller’s  Box  scheme  [37]  and  is  applied 
to  the  asymptotic  form  of  the  Navier-Stokes  equations 
(the  sublayer  equations)  in  a  thin  region  (typically  less 


than  1%  of  the  boundary  layer)  adjacent  to  solid  sur¬ 
faces.  The  sublayer  approximation  has  been  validated 
for  both  2-D  [38]  and  3-D  [30]  applications.  The  ex¬ 
plicit  algorithm  of  MacCormack  [39]  is  applied  to  the 
full  Navier-Stokes  equations  on  grid  points  exterior  to 
the  computational  sublayer.  This  includes  nearly  all  of 
the  boundary  layer  and  the  external  inviscid  re^on.  The 
code  has  been  employed  successfully  to  predict  turbulent 
interactions  fbr  the  3-D  sharp  fin  [14],  swept  compres¬ 
sion  comer  [25,  26]  and  the  double  fin  configuration  [30]. 
Computations  described  in  this  paper  were  performed 
either  on  the  now  decommissioned  4-pipe  CYBER  205 
at  the  NASA  Ames  Research  Center  or  on  the  2-pipe 
VPS-32  machine  at  NASA  Langley  Research  Center. 

3  Classification  of  Computations 

Two  configurations  are  utilised  in  this  research  — 
the  SFC  (Fig.  1)  and  the  SDFC  (Fig.  4).The  single-fin 
configuration  (SFC)  is  employed  to  study  the  effects  of 
auction  on  the  flowfield.  None  of  the  computations  on 
the  symmetric  double-fin  configuration  (SDFC)  incor¬ 
porate  suction.  The  attempt  at  control  arises  rather 
&om  the  geometry  itself  since  the  vorticity  generated 
by  the  two  fins  tends  to  counteract.  For  the  purposes  of 
the  following  discussion,  the  x-axis  points  in  the  general 
streamwise  direction,  yia  plate  normal  and  z  is  the  span- 
wise  axis.  Also,  the  theoretical  boundary  layer  thick¬ 
ness  is  utilised  to  scale  distances  for  Cases  A,  B,  C,  F 
and  G  (see  below).  The  experimental  value  is  utilized 
for  Cases  D  and  E  (Table  1). 

3.1  Suction  on  SFC 

In  each  case,  a  particular  area  (see  below)  is  de- 
mazkated  and  the  mass  flux  at  each  grid  point  in  the 
area  is  specified  as  a  ratio  of  the  upstream  freestream 
mass  flux.  The  bleed  parameter,  B.P.  I,defined  as  fol¬ 
lows: 

Bleed  Parameter  (B.P.)  I  =  ^  *  (1) 

PooUao 

is  utilised  to  characterise  the  magnitude  of  suction.  In 
this  equation,  p  is  the  density,  w  and  v  are  the  Cartesian 
components  of  the  velocity  in  the  streamwise  and  plate 
normal  directions  and  subscripts  in  and  oo  indicate  wall 
and  fieestream  respectively.  The  implementation  corre¬ 
sponds  to  porous  bleed  as  opposed  to  slot  bleed.  The 
parameter  B.P.  I  does  not  incorporate  the  area  of  bleed 
and  therefore,  provides  only  a  local  estimate  of  bleed 
magnitude.  A  more  physical  value  may  be  computed  by 
approximating  the  percentage  of  fluid  in  the  incoming 
boundary  layer  bled  off.  This  requires  an  estimate  of  the 
amount  of  incoming  flow  that  actually  “sees”  the  bleed 
region.  This  guess  may  be  made  in  a  simple  manner 
by  simply  projecting  the  ray  indicating  start  of  suction 
on  the  upstream  end  of  the  domain.  For  example,  for 
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Legend: 


6  -  Boundary  layer  thickness 
IL,JL,KL  -  Points  in  a,y,s  directions 
NSLi  -  Grid  points  in  fin  sublayer 
Re  -  Reynolds  number 
SFC  -  Single  fin  configuration 
SYS  n  -  2-Pipe  VPS-32  (Langley) 

T|  -  Total  temperature 


B.P.  -  Bleed  Parameter 

M  -  Mach  number 

NSLa  -  Grid  points  in  plate  sublayer 

Pt  -  Total  pressure 

SYS  I  -  4-Pipe  CYBER  205  (Ames) 

Tc  -  Characteristic  time  (=  L/Uoo) 

Tw  -  Wall  temperature  (fin  and  plate) 


Table  1:  Classification  of  Computations 


Requirement 

*1+  <3 

*1+R<3 

NPBL 

Avg. 

■^1 

Avg. 

i  Avg. 

Max 

Avg. 

Max 

Case  A 

1.47 

3.03 

1.47 

3.03 

47.4 

97.4 

16 

Case  B 

1.47 

3.01 

1.46 

3.01 

I  0.30 

0.41 

50.4 

97.0 

16 

Case  C 

1.99 

4.21 

1.62 

3.03 

SOU 

KS3II 

59.3 

135.3 

16 

Case  D 

0.89 

2.10 

0.89 

2.10 

0.0 

0.0 

28.6 

69.4 

16 

Case  £ 

1.10 

5.02 

0.91 

2.18 

1.42 

1.44 

35.26 

161.6 

16 

Case  F 

0.91 

1.01 

0.91 

1.01 

0.0 

KSH 

19.1 

21.2 

16 

Case  G 

1.05 

1.20 

1.05 

1.20 

HESI 

KSfll 

19.4 

24.3 

16 

Legend:  -  First  grid  point  from  boundary  in  wall  units 

N  -  Bleed  correction  factor  of  Cebeci 
m  -  Local  bleed  magnitude 
/iw  -  Molecular  viscosity  at  wall 

-  Height  of  sublayer  in  wall  units 
NPBL  -  Points  in  boundary  layer 


Table  2:  Grid  resolution  requirements 
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Region  I  in  Fig.  3,  the  boundary  layer  in  distance  Li 
may  be  assumed  to  “see”  the  bleed  region.  B.P.  11  may 
be  expressed  as: 

B.p.n=.  (2) 

^  ~  S*)p„u^L 

where  6  is  the  boundary  layer  thickness,  6*  is  the  dis¬ 
placement  thickness,  At  is  the  area  of  bleed  and  L  is 
the  appropriate  upstream  distance.  Two  fin  angles  are 
utilised  for  bleed  computations: 

20  deg  SFC 

Case  A  corresponds  to  flow  past  the  20  SFC  with¬ 
out  bleed.  For  Cases  B  and  C,  1%  (=  B.P.  I)  is  imposed 
on  Redons  I  and  U  (Fig.  3)  respectively.  The  first  re¬ 
gion  (Region  I  in  Fig.  3)  extends  between  the  line  of 
upstream  influence  and  the  TISF.  The  second  region 
(Region  II)extends  between  the  TISF  and  the  fin  base. 
The  bleed  and  flow  parameters  are  presented  in  Table  1. 
Note  that  Cases  A  and  B  have  already  been  described 
in  [30]  (identified  there  as  the  “no  bleed”  and  “1%  bleed” 
cases)  and  are  included  here  only  for  comparison  with 
Case  C. 

8  deg  SFC 

The  recent  experimental  investigation  by  Barn¬ 
hart  et  al  [31]  on  the  effect  of  suction  on  supersonic 
flow  past  an  8”  SFC  is  numerically  investigated  in  this 
research.  The  flow  parameters  are  summarised  in  Tk- 
ble  1.  Two  computations  are  performed  on  this  geom¬ 
etry,  a  no  bleed  case  (Case  D)  and  a  case  incorporating 
suction  on  the  area  displayed  in  Fig.  5  (Case  E).  As  for 
the  20“  SFC,  bleed  is  imposed  by  specifying  the  bleed 
mass  flux  on  the  plate  at  each  grid  point  in  the  bleed 
region.  The  relevant  bleed  parameters  are  provided  in 
Ikble  1. 

3.2  Flow  past  the  SDFC 

Computations  are  performed  on  two  distinct  SD- 
FCs  (Fig.  4)  with  fin  angles  4x4  (Case  F)  and  8x8 
(Case  G).  The  overall  flow  parameters  are  listed  in  Ta¬ 
ble  1.  The  SDFC  has  two  geometrical  length  scales  — 
the  sepuation  between  the  fins  (£«  in  Fig.  4)  and  the 
length  of  the  fins  (X/tn)*  L^  is  chosen  to  be  33.33^  in 
accordance  with  experimental  considerations.  The  value 
of  L fin  (696)  is  chosen  such  that  for  Case  F,  the  inviscid 
secondary  shock  hits  the  fin  at  the  trailing  edge  of  the 
domain. 

4  Boundary  Conditioru 

The  boundary  conditions  employed  for  both  con¬ 
figurations  are  described  together  with  reference  to 
Figs.  1  and  4.  The  precise  mapping  from  physical  to 


computational  domain  is  achieved  with  a  computational 
grid  described  briefly  in  Section  5  and  in  detail  in  [40]. 

•  Upstream  boundary  [ABHG]:  The  flow  at  this 
boundary  (located  56  upstream  of  the  fin  leading 
edge  for  Cases  A,B,C,F  &  G,  4.66  for  Cases  D 
b  E)is  assumed  to  be  two-dimensional  and  is  gen¬ 
erated  with  a  separate  2-D  compressible  boundary 
layer  code  [36]  utilising  the  Baldwin-Lomax  turbu¬ 
lent  eddy  viscosity  model.  For  the  SFC  computa¬ 
tions,  the  momentum  thickness  is  matched  while  for 
the  SDFC  computations,  the  anticipated  boundary 
layer  thickness  is  matched. 

•  Plane  of  symmetry  upstream  of  fin  leading  edge 
[AFLG]:  S3rmmetry  boundary  conditions  are  ap¬ 
plied  at  this  boundary  i.e.,  the  normal  component 
of  the  velocity  and  the  normal  derivatives  of  the  re¬ 
maining  flow  quantities  are  set  to  sero.  For  both 
the  SFC  and  the  SDFC, 


«,  =  0; 

dh 


=  0 


(3) 


where,  p  is  the  density,  v,  «  and  w  are  the  x,y,  and 
z  components  of  the  velocity  vector  (=  {u, «,  w}), 
n  is  the  normal  to  the  boundary,  p  and  T  are  the 
pressure  and  temperature  respectively. 

•  Fin  surface  [LFKE]:  Since  this  is  a  solid  surface,  the 
velocity  vector  and  the  normal  pressure  gradient 
are  assumed  sero  and  a  fixed  surface  temperature 
is  specified  i.e.. 


pv  =  Qi  T  =  T*;  ^  =  0  (4) 

where  the  subscript  it;  refers  to  wall  conditions.  The 
thermal  boundary  conditions  employed  for  all  cases 
are  presented  in  Table  1. 

•  Flat  plate  [ABCDEF]:  The  two  components  of  the 
velocity  vector  along  the  flat  plate  are  specified 
sero.  Bleed  is  applied  in  the  appropriate  areas 
for  Cases  B,  C  and  E  by  specifying  the  quantity 
PWVw/pooVoo  (B.P.  I  in  Table  1).  The  normal  pres¬ 
sure  gradient  is  assumed  sero  and  the  surface  tem¬ 
perature  is  fixed  as  for  the  fin  surface  (Table  1). 


pv  =  m;  u 


ii;  =  0;  r=T,; 


8n 


0 


(5) 


•  Downstream  boundary  [EKJD]:  In  every  case,  the 
streamwise  length  of  the  physical  domain  is  re¬ 
stricted  so  that  the  outgoing  flow  is  supersonic  (ex¬ 
cept  close  to  the  flat  plate,  of  course)  and  as  such, 
the  sero  gradient  extrapolation  condition  {d/8z  = 
0)  is  applied.  For  Cases  A,B  and  C,  this  boundary 
is  located  266  downstream  of  the  fin  leading  edge. 
For  Cases  D  and  E  this  distance  is  17.756  and  for 
Cases  F  and  G,  the  value  is  696. 
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•  Right  boundaiy  [BCDJIH]:  This  boundsty  condi¬ 
tion  is  treated  separately  for  the  SFC  and  the 
SDFC: 


-  SFC:  The  boundary  is  assumed  to  be  suffi¬ 
ciently  far  from  the  fin  to  ensure  that  the 
boundary  layer  is  locally  two-dimensional  and 
therefore  a  simple  gradient  boundary  condition 
djdz  =  0  is  employed.  This  requires  that  the 
shock  wave  pass  out  the  downstream  boundary 
[EKJD]. 

-  SDFC:  This  is  a  symmetry  boundary  and  the 
same  conditions  are  applied  as  at  the  plane 
of  symmetry  upstream  of  the  fin  leading  edge 
[AFLG]. 


•  Top  Boundary  [HUKLG]:  This  boundary  is  lo¬ 
cated  far  enough  (roughly  8  j)  away  from  the  corner 
(formed  by  the  plate  and  the  fin)  to  permit  the  as¬ 
sumption  of  two-dimensional  fiow  in  planes  parallel 
to  the  boundary.  The  boundary  conditions  applied 
are  therefore: 

«» 


All  boundary  derivatives  are  implemented  in  single¬ 
sided  first  order  accurate  form. 


S  Computational  Details 

Non-uniform  Cartesian  grids  are  employed  in  all 
calculations.  The  3-D  grid  is  composed  of  a  set  of 
streamwise  (z  =  constant)  planes  (equally  spaced  with 
An  =  ^  for  all  cases  except  D  and  E).  Within  each 
plane,  two  separate  grid  systems  are  employed  —  the 
"ordinary”  and  the  "sublayer”  grid  [33].  The  ordinary 
grid  points  a4j&cent  to  the  plate  and  the  fin  are  dis¬ 
tributed  using  a  combination  of  geometrically-stretched 
and  uniform  spacings  in  the  general  fashion  of  [20].  The 
number  of  grid  points  utilised,  CPU  and  corresponding 
fiow  development  time  required  for  each  case  are  pre¬ 
sented  in  Table  1. 

Adequate  resolution  of  different  regions  of  the  fiow 
impose  various  constraints  on  permissible  grid  sises.  A 
complete  set  of  criteria  forjudging  the  acceptability  of  a 
grid  system  do  not  exist  at  present,  except  perhaps  for 
the  process  of  successive  grid  refinement  which  is  gener¬ 
ally  not  feasible  for  3-D  Navier  Stokes  calculations.  The 
determination  of  the  acceptability  of  a  given  grid  sys¬ 
tem,  must  be  based  on  experience  and  previous  numer¬ 
ical  computations  [33].  The  characteristics  of  the  grids 
employed  for  every  case  satisfy  constraints  described  in 
the  literature  in  the  average  sense  (see  Table  2).  The 
maximum  values  are  localised  in  a  few  locations  such  as 
in  bleed  regions  and  at  the  fin  leading  edge,  where  the  lo¬ 
cal  skin  friction  r«  shows  steep  increases  (see  Fig.  11  for 
example),  making  it  unfeasible  to  satisfy  all  constraints 
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5.81 

14 

1.13 
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Table  3:  Location  of  Experimental  Surveys 


locally.  Previous  computations  with  similaz  values  have 
compared  well  with  experiment  [33]. 

6  Results 

For  the  purpose  of  brevity,  several  results  are  pre¬ 
sented  without  figures.  A  more  detailed  description  of 
this  work  may  be  found  in  [40]. 

6.1  Suction  on  20  deg.  SFC 

Computed  results  are  first  formally  compared  with 
available  experimental  data  in  the  form  of  boundary 
layer  pitot  pressure  and  yaw  angle  profiles  and  surface 
pressures.  Since  no  bleed  is  imposed  in  the  experiments, 
direct  comparison  with  theory  is  limited  to  Case  A  (no 
bleed).  In  the  following  discussion,  the  upstream  bound¬ 
ary  layer  thickness  £ge  (theoretical)  is  utilised  to  scale 
distances  (c£  [13]  for  alternate  scales).  Streamwise  dis¬ 
tances  (a),  unless  otherwise  mentioned,  are  measured 
relative  to  the  location  of  the  inviscid  shock  (a«)  at  the 
relevant  spanwise  (x)  position.  The  distance  normal  to 
the  flat  plate  is  denoted  y  .  Flow  variables  are  normal¬ 
ised  by  their  freestream  values  upstream  of  the  interac¬ 
tion. 

Effect  of  suction  on  Flow  Variables 

The  fourteen  experimental  stations  of  Shapey  and 
Bogdonoff  [14]  are  located  in  the  general  shape  of  a  cross 
(Fig.  6)  with  stations  1  and  2  located  upstream  of  the 
interaction,  3  and  13  at  about  the  line  of  upstream  in¬ 
fluence,  stations  4,  5,  6  and  12  lie  in  Region  I,  while 

7  through  11  and  14  lie  in  Region  II.  The  coordinates 
of  these  stations  are  provided  in  Table  3.  Experimen¬ 
tal  and  numerical  results  available  in  the  form  of  pitot 
pressure  and  yaw  angle  boundary  layer  surveys  are  com¬ 
pared  at  each  of  the  above  stations,  although,  for  the 
purpose  of  brevity,  only  some  plots  (stations  3,  S,  and  8) 
are  presented  here. 

At  stations  1  and  2  located  upstream  of  the  in¬ 
teraction,  the  pitot  pressure  and  yaw  angle  profiles  (not 
shown  —  see  [40])  are  unaffected  by  suction  in  either  re¬ 
gion  and  resemble  profiles  for  2-D  boundary  layer  flows 


6 


aa  expected.  The  pitot  pieMuie  and  yaw  angles  at  sta¬ 
tion  3  are  shown  in  Fig.  7  wheie  the  ordinate  {y/6)  is  the 
distance  normal  to  the  flat  plate  and  the  abscissa  is  the 
pitot  pressure  (Fig.  7  (a)),  normalised  by  the  freestream 
pitot  pressure,  or  yaw  angle  in  degrees  (Fig.  7  (b)).  A 
perceptible  drop  in  pitot  pressure  because  of  bleed  is  ev¬ 
ident  at  about  y/j  =  0.1  for  Region  I  bleed  in  contrast 
to  Region  II  bleed  which  shows  negligible  influence  at 
this  station.  A  similar  conclusion  may  be  drawn  from 
the  yaw  angle  profile.  Region  I  bleed  displays  reduced 
yaw  angles  up  to  about  y/^oo  =  0.8  while  Region  n 
bleed  has  no  significant  effect.  This  general  behavior 
is  also  seen  at  stations  4  (not  shown)  and  5  (Fig.  8) 
both  of  which  lie  upstream  of  Region  IL  At  both  these 
stations,  whereas  Region  II  bleed  has  no  infinence.  Re¬ 
gion  I  bleed  exhibits  a  "cross-over**  point,  i.e.,  near  the 
flat  plate  (up  to  about  y/Sao  =  0.7  for  station  5)  the 
pitot  pressure  increases  above  its  no  bleed  value  while, 
at  larger  distances  from  the  plate,  the  pitot  pressure 
is  slightly  lower.  The  yaw  angles  at  both  stations  (4 
and  5)  are  lower  for  Region  I  bleed  with  no  effect  of  Re¬ 
gion  II  bleed.  At  station  8  (Fig.  9),  the  pitot  pressure  for 
Region  II  bleed  is  larger  very  close  to  the  flat  plate  ac¬ 
centuating  the  previously  observed  t  shaped  structure. 
In  the  major  portion  of  the  boundary  layer  however,  the 
pitot  pressure  closely  parallels  the  no  bleed  profile.  Yaw 
angles  for  Region  II  bleed  show  no  significant  deviation 
from  no  bleed  (Case  A)  values  except  in  the  free  stream. 

The  computed  wall  pressure  is  compared  with  the 
experimental  data  of  Goodwin  [10]  in  Fig.  10  at  the 
spanwise  location  z/6  =  5.  The  ordinate  is  the  static 
pressure  normalised  by  the  freestream  static  pressure. 
The  abscissa  is  the  streamwise  distance  (measured  rela¬ 
tive  to  the  inviscid  shock)  normalised  with  6^  .  Com¬ 
parison  of  Case  A  with  experimental  data  indicates  good 
agreement  with  regard  to  upstream  infiuence.  The  ex¬ 
perimental  data  indicates  a  plateau  region  at  about 
—  0)  which  is  not  adequately  reproduced  in  the 

computation  perhaps  for  reasons  of  grid  resolution.  The 
total  pressure  rise  equals  the  expected  theoretical  value 
in  the  invisced  region  and  is  independent  of  bleed  region 
and  magnitude.  The  effect  of  Region  I  bleed  (Case  B) 
is  evidently  to  reduce  upstream  influence  thus  reducing 
modestly  the  distance  over  which  the  pressure  rise  is 
achieved.  Region  II  bleed  (Case  C)  does  not  influence 
the  upstream  influence  or  the  total  pressure  rise. 

The  localised  effect  of  bleed  on  skin  friction  co¬ 
efficient  is  exhibited  in  Fig.  11  at  x/6  =  5.  For  both 
bleed  regions,  local  skin  friction  values  display  steep  in¬ 
creases.  It  is  apparent  that  a  significant  drag  penalty  is 
associated  with  the  employment  of  suction. 

Fig.  12  displays  static  pressure  contours  on  the 
z/6  =  5.81  plane.  The  high  concentration  of  static 
pressure  contours  in  the  region  away  from  the  plate 


(y/j  >  3)  indicate  the  existence  of  the  inviscid  shock 
which  is  smeared  over  roughly  three  (3)  grid  points  at 
this  section.  Near  the  plate,  the  shock  wave  forms  two 
compression  fans  -  an  upstream  (left)  and  a  downstream 
(right)  fan.  The  effect  of  bleed  in  Region  I  is  evidently 
to  reduce  the  sise  of  the  upstream  fan.  Region  II  bleed, 
on  the  other  hand,  has  no  significant  effect  on  either  fan. 

Effect  of  Sneiion  on  Flowiield  Stractnre 

The  computed  solutions  are  utilised  to  examine 
the  effect  of  bleed  on  the  flow  structure.  Fig.  13  dis¬ 
plays  surface  skin  friction  lines  (often  denoted  "surface 
streamlines")  for  Cases  A,  B  and  C  respectively.  These 
are  generated  with  a  particle  tracing  code  with  the  ap¬ 
propriate  surface  stresses  replacing  velocities.  The  loca¬ 
tion  of  the  inviscid  shock  is  indicated  in  each  figure.  The 
line  of  coalescence  clearly  persists  for  all  cases.  Whereas 
Region  I  suction  (Fig.  13(b))  tends  to  align  the  line  of 
coalescence  with  the  theoretical  inviscid  line,  Re^on  II 
suction  (Fig.  13(b))  has  no  significant  influence  despite 
the  previously  discussed  steep  local  increase  in  skin  fric¬ 
tion  values.  This  is  another  manifestation  of  the  con- 
clurion  derived  from  the  surface  pressure  plots,  vis..  Re¬ 
gion  II  bleed  has  no  effect  on  upstream  influence. 

A  further  analysis  of  the  computed  flowfield  is  car¬ 
ried  out  with  particle  traces  employing  a  3-D  code  [41]. 
Particles  are  "released”  at  various  heights  above  the  flat 
plate  at  the  upstream  boundary  and  their  paths  are  plot¬ 
ted  to  study  gross  flow  features.  Fig.  14  shows  traces 
of  particles  released  at  y/S^  —  0.0001,  within  the  sub¬ 
layer.  For  Case  A,  these  particles  coalesce  to  form  the 
line  of  coalescence  observed  in  the  surface  streamlines. 
For  Case  B,  the  particles  are  ingested  by  the  plate  in  the 
region  of  bleed.  Particles  for  Case  C  follow  trajectories 
similar  to  those  for  Case  A.  Fig.  15  shows  traces  of  par¬ 
ticles  released  at  y/^oo  =  •  For  both  Cases  A  and  B, 

the  particles  are  entrained  in  the  vortical  structure.  For 
Case  C,  those  particles  originating  near  the  leading  edge 
of  the  fin  at  the  upstream  boundary  (small  values  of  the 
spanwise  coordinate  z),  are  ingested  by  the  fiat  plate  in 
Region  II  (note  the  trajectories  of  the  first  seven  par¬ 
ticles  from  the  fin  leading  edge).  The  other  particles 
are  swept  into  the  vortical  structure  as  for  the  other 
cases.  It  is  evident  that  only  those  particles  whose  no 
bleed  trajectory  passes  close  to  the  flat  plate  (to  within 
a  certain  height  whose  value  is  difficult  to  estimate)  are 
ingested.  To  study  the  origin  of  the  flow  that  is  ingested 
by  the  plate  in  Case  5,  a  series  of  particles  were  released 
near  the  plate  (in  the  bleed  region)  and  their  motion 
was  traced  in  reverse.  The  traces  are  shown  in  Fig.  16. 
These  particles  generally  originate  near  the  fin  leading 
edge.  In  the  no  bleed  case,  they  would  then  have  then 
risen  up  and  around  the  vortical  core.  The  above  figures 
show  clearly  that  the  mean  flowfield  pattern  continues 


7 


to  be  dominated  by  the  vortical  structure  for  both  bleed 
re^ons. 

6.2  Suction  on  8  deg  SFC 

The  geometrical  configuration,  bleed  region  and 
other  significant  parameters  of  the  two  computations 
(Cases  D  and  E)  on  the  8"  SFC  are  described  in  Fig.  5 
and  Table  1.  As  for  the  20*  SFC,  results  ate  first  com¬ 
pared  with  experimental  data  and  the  flow  is  then  ex¬ 
amined  with  particle  traces. 

Compimison  with  Experiment 

The  experimental  data  utilised  for  comparison 
purposes  was  recorded  at  the  NASA  Lewis  Research 
Center  and  is  repotted  in  Barnhart  et  al  [31].  They 
utilise  two  Mach  numbers,  2.48  and  2.95.  For  the  lat¬ 
ter,  relevant  to  this  research,  results  are  published  only 
in  the  presence  of  bleed  and  include  surfSsce  (static  pres¬ 
sures  and  wall  flow  angles)  and  boundary  layer  measure¬ 
ments  (pitot  pressure  and  yaw  angles).  Each  of  these  is 
compared  below  with  computed  results. 

It  is  important  to  note  one  fundamental  distinc¬ 
tion  between  the  computed  and  experimental  configu¬ 
ration.  In  the  experiment,  suction  is  applied  by  per¬ 
forations  on  the  plate  in  coiyunction  with  a  vacuum 
boost  pump  downstream  of  the  suction  plenum.  Since 
the  surface  pressure  increases  into  the  interaction,  the 
experimental  bleed  rate  is  variable,  increasing  generally 
in  the  streamwise  direction  at  any  given  spanwise  lo¬ 
cation.  For  lack  of  precise  information  on  local  bleed 
variation  on  the  plate,  the  computation  assumes  a  uni¬ 
form  bleed  mass  rate  such  that  the  total  bleed  mass 
rate  over  the  entire  bleed  area  matches  experiment.  It 
is  expected  therefore,  that  computed  bleed  values  ex¬ 
ceed  experimental  values  in  the  vicinity  of  the  upstream 
end  of  the  bleed  domain.  The  reverse  situation  may  be 
expected  at  the  downstream  end.  In  the  following  dis¬ 
cussion,  distances  are  normalised  with  the  ezpenmen- 
toi  upstream  boundary  layer  thickness  (S  =  3.112cm.). 
The  term  “experimental  data”  is  assumed  to  imply  in 
the  pretence  of  bleed  since  no  data  is  available  without 
suction  for  this  configuration.  The  strict  comparison  is 
therefore  restricted  to  Case  E. 

The  computed  (Case  E  only)  and  experimental 
surface  pressures  at  four  spanwise  locations  are  plot¬ 
ted  in  Fig.  17.  Note  that  the  profiles  at  all  locations 
collapse  roughly  on  a  single  curve.  Further,  while  the 
upstream  pressures  ((a  —  x,)/6  ~  —4)  match  well  as 
expected,  the  computation  displays  significantly  lower 
upstream  influence.  Experimental  pressures  begin  to 
rise  at  (a  —  a,)/6  ~  -1.5  whereas  the  initial  com¬ 
puted  pressure  rise  occurs  at  {z  -  x,)/S  ~  0,  at  the 
theoretical  inviseid  shock  location.  This  discrepancy  is 
partly  attributable  to  the  assumption  of  uniform  bleed 


in  the  computation.  Since  the  effect  of  bleed  is  to 
lower  upstream  influence,  a  higher  local  computed  bleed 
value  (over  experiment)  may  be  expected  to  result  in 
a  lower  upstream  influence.  The  maximum  deviation 
25%)  occurs  near  the  location  of  the  inviseid  shock 
((a  —  Zt)/S  ~  0)  and  beyond.  Near  the  center  of  the 
suction  region,  the  computed  and  experimental  values 
match,  resulting  from  a  higher  rate  of  pressure  rise  for 
the  former.  The  overall  pressure  rise  at  the  spanwise 
location  (Z  =  2.46)  matches  experiment  and  is  only 
slightly  lower  (by  ^  3%)  than  the  theoretical  pressure 
ri«e  (Pj/Pi  =  1*78). 

The  comparison  with  experimental  surface  yaw 
values  at  the  same  spanwise  locations  is  displayed  in 
Fig.  18.  The  value  ct/6  is  plotted  on  the  ordinate  where 
a  is  the  turning  angle  and  9  is  the  fin  angle.  Upstream 
of  the  interaction  ((z  —  z,)/6  —4),  yaw  angles  are 

roughly  sero.  l^thin  the  interaction  (—3  <  {z—x,)IS  < 
2),  experimental  values  display  considerable  scatter  and 
indicate  a  more  gradual  turning  of  the  flow  than  does 
the  computation.  The  discrepancy  is  again  partly  at¬ 
tributable  to  differences  in  local  bleed  rates.  Computa¬ 
tions  predict  lower  turning  angles  at  the  upstream  edge 
of  the  bleed  region  because  local  bleed  values  are  higher 
than  probably  exist  in  experiment.  Both  computation 
and  experiment  display  maximum  turning  angles  lower 
than  the  shock  angle  (corresponding  to  the  value  3.23  in 
Fig  18)  although  maximum  experimental  values  exceed 
computed  by  roughly  15%.  Note  that  maximum  turning 
occurs  behind  the  inviseid  shock.  Further  downstream, 
computed  values  display  turning  angles  twice  the  fin  an¬ 
gle  whereas  experimental  values  approach  the  fin  angle. 

Fig.  19  exhibits  the  location  of  the  nine  experi¬ 
mental  boundary  layer  surveys  where  pitot  pressure  and 
yaw  angle  profiles  were  measured.  Fig  20  displays  pitot 
pressure  profiles  at  each  of  the  nine  stations.  At  sta¬ 
tion  1,  upstream  of  the  interaction,  the  experimental 
values  (markers)  match  both  the  no  bleed  (solid  line) 
and  bleed  (dotted  line)  cases  and  resemble  a  2-D  pro¬ 
file  similar  to  that  existing  at  the  upstream  end  of  the 
domain.  At  station  2,  the  curves  are  generally  similar 
to  those  at  station  1  although  in  the  region  between 
0.3  <  y/6  <  1,  experimental  values  are  slightly  lower. 
At  the  upstream  edge  of  the  bleed  domain  (station  3), 
both  computed  and  experimental  values  in  the  presence 
bleed  display  modestly  lower  values  in  the  upper  half 
of  the  boundary  layer.  At  station  4,  on  the  shock  loca¬ 
tion,  Case  D  displays  a  small  pitot  pressure  overshoot. 
The  effect  of  bleed  is  to  increase  pitot  pressure  values 
up  to  about  y/6  ~  0.3  compared  to  the  no  bleed  case, 
beyond  which  the  values  are  lower.  A  simUar  observa¬ 
tion  was  made  previously  for  the  20*  SFC  (see  Fig.  8  for 
example).  The  experimental  results  agree  with  Case  E 
beyond  y/6  >  0.3  although  closer  to  the  plate,  more 
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confotmity  is  obt&ined  with  the  no  bleed  case.  A  simi¬ 
lar  crossover  of  no  bleed  and  bleed  profiles  is  observed 
at  station  5,  the  first  station  located  downstream  of  the 
theoretical  shock.  Comparing  Case  E  with  experiment, 
in  the  region  Q  <  y/6  <  0.8,  the  computed  values  over¬ 
predict  pitot  pressure  values.  The  situation  is  reversed 
higher  into  the  boundary  layer.  At  stations  6,  T,  8  and 
9,  the  profiles  all  display  similar  qualitative  behaviour. 
The  “crossover”  of  bleed  and  no  bleed  profiles  occurs 
successively  higher  into  the  boundary  layer.  In  gen¬ 
eral,  experimental  values  agree  with  Case  E  close  to  the 
plate  although  higher  into  the  boundary  layer,  experi¬ 
mental  values  match  no  bleed  results  quite  closely,  indi¬ 
cating  that  the  effect  of  bleed  propagates  higher  into  the 
boundary  layer  in  the  computation  than  in  experiment. 

The  computed  and  experimental  yaw  angles  are 
disifiayed  in  Fig.  21.  At  station  1,  upstream  of  the 
interaction,  yaw  angles  are  negligible  throughout  the 
boundary  layer  as  anticipated.  At  station  2,  the  effect 
of  bleed  is  to  lower  yaw  angle  values  below  the  no  bleed 
values  up  to  about  y/S  ^  0.5.  At  the  surface,  computa¬ 
tional  values  exceed  experimental  values.  At  the  start 
of  the  bleed  region  (station  3),  computed  yaw  angles 
for  Case  E  assume  small  values  through  out  the  bound¬ 
ary  layer  although  in  the  absence  of  bleed  (Case  D), 
the  flow  exhibits  significant  turning  at  this  station.  Ex¬ 
perimental  values  lie  between  those  for  Cases  D  and  E. 
Downstream  of  the  inviscid  shock  (station  5),  the  cor¬ 
respondence  between  Case  E  and  experiment  is  good 
with  the  computations  tending  to  slightly  underpredict 
yaw  angles.  The  effect  of  bleed  is  to  reduce  yaw  angles 
in  the  boundary  layer  by  as  much  as  ~  50%  at  the  sur¬ 
face.  Further  downstream,  at  station  6,  although  Case  E 
displays  lower  yaw  angles,  experimental  values  paral¬ 
lel  Case  D  more  closely.  At  stations  7  and  8,  located 
downstream  of  the  bleed  region,  the  effect  of  suction  is 
more  modest.  The  comparison  between  experiment  and 
Case  E  is  very  good  at  station  7.  Experimental  obser¬ 
vations  display  lower  turning  at  close  proximity  to  the 
plate  at  station  8.  At  station  9,  located  far  downstream 
of  the  bleed  region,  the  effect  of  bleed  is  restricted  to 
the  region  0  <  y/5  <  0.6.  The  experimental  values, 
which  display  significantly  lower  values  at  this  location, 
are  suspect  since  the  angle  obtained  in  the  free  stream 
is  far  lower  than  anticipated  (equal  to  the  fin  angle  since 
this  station  is  downstream  of  the  shock). 

Analysis  of  Computed  Flowfield 

The  overall  effects  of  bleed  are  observed  to  be  the 
same  as  for  the  previously  discussed  20*  SFC.  Fig.  22 
compares  surface  static  pressure  at  two  spanwise  loca¬ 
tions  for  Cases  D  and  E.  Upstream  influence,  determined 
by  initial  pressure  rise,  is  reduced  and  the  rate  of  pres¬ 
sure  rise  in  the  interaction  increases  with  bleed.  The 


total  pressure  rise  remains  the  same.  Local  skin  &iction 
values  (not  shown)  display  steep  increases  in  the  bleed 
repon  and  are  not  affected  elsewhere.  The  effect  of  bleed 
on  the  shock  structure  is  also  similar  to  that  observed 
for  the  20*  fin  with  Region  I  bleed.  The  effect  of  bleed 
is  to  completely  eliminate  the  upstream  fan  resulting 
in  a  shock  structure  that  appears  entirely  inviscid  (not 
shown).  No  significant  influence  of  bleed  is  observed  in 
the  inviscid  regions. 

The  computed  surface  streamlines  are  displayed 
in  Figs.  23  and  24  for  Cases  D  and  E  respectively.  For 
Case  D,  surface  streamlines  display  maximum  turning 
angles  roughly  equal  to  the  shock  angle.  No  clear  line 
of  coalesence  is  visible  though  the  lines  display  bunch¬ 
ing  slightly  upstream  of  the  inviscid  shock  trace.  For 
Case  E,  significant  reduction  in  surface  angularity  is  ob¬ 
served  in  the  bleed  region  only.  Streamlines  originating 
near  the  fin  leading  edge  follow  trajectories  similar  to 
those  for  Case  D.  Upon  reaching  the  bleed  repon  how¬ 
ever,  they  straighten  out  and  move  parallel  to  the  up¬ 
stream  flow  for  a  short  distance  gradually  obtaining  the 
shock  angle  downstream  of  the  bleed  region.  Several 
lines  passing  through  the  furthermost  upstream  end  of 
the  bleed  region  form  a  single  group  of  closely  bunched 
lines  passing  out  the  extreme  downstream  end  of  the 
bleed  region. 

The  flowfield  is  analysed  with  exhaustive  sets  of 
particle  traces  released  at  various  heights  above  the 
plate.  For  brevity,  only  a  few  sets  are  discussed.  Parti¬ 
cles  released  within  the  sublayer  {jy/S  =  0.0001)  follow 
the  general  trajectory  of  the  surfhce  streamlines  as  an¬ 
ticipated  (see  Fig.  25  (a))  for  Case  D.  In  the  presence 
of  bleed  however,  particles  are  ingested  (Fig.  25  (b))  at 
the  leading  edge  of  the  bleed  region  if  released  appro¬ 
priately  close  to  the  fin.  Particles  at  larger  spanwise 
distances  escape  the  influence  of  suction.  A  similar  pat¬ 
tern  is  observed  up  to  about  y/S  =  0.01.  Further  into 
the  boundary  layer,  at  y/S  =  0.1  (Fig.  26),  Case  D  dis¬ 
plays  significant  flow  turning  although  the  clear  vortical 
structure  observed  for  the  20*  SFC  does  not  exist.  Sev¬ 
eral  particles  originating  near  the  fin  and  whose  no  bleed 
trajectory  passes  over  the  bleed  region  are  ingested  in 
Case  E  roughly  at  the  center  of  the  bleed  region.  Par¬ 
ticles  released  further  from  the  plate  (y/5  =  0.2)  are 
ingested  at  the  trailing  edge  of  the  bleed  region  (not 
shown).  Particles  released  at  larger  spanwise  distances 
are  not  affected.  At  y/S  —  0.3  (Fig.  27),  none  of  the 
particles  are  ingested  although  the  traces  are  slightly 
altered  by  suction  and  display  modestly  lower  yawing. 
At  y/5  =  1  (Fig.  28),  the  effect  of  bleed  is  negligible. 

6.3  Flow  past  double  fin  configurations 

Both  computations  on  the  SDFC  (see  Table  1)  are 
described  together.  Since  experimental  data  is  unavail- 
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able  at  thia  time,  only  computed  leaults  ate  diacuaaed. 
Note  that  the  domain  length  ia  identical  for  both  the 
weaker  (4  x  4)  and  the  ationger  (8  x  8)  inteiactiona. 
Hence,  the  free  atieam  flow  paaaea  through  a  larger  num¬ 
ber  of  ahock  wavea  for  the  latter  caae. 

Fig.  29  shows  computed  surface  pressure  contours 
plotted  after  symmetrically  extending  the  flowfield.  The 
theoretical  pressure  rise  after  each  shock  for  both  cases 
is  presented  in  Table  4.  These  figures  are  qualitatively 
similar  for  the  initial  (primary)  pressure  rise  though  the 
pressure  rises  faster  for  the  larger  angle  as  is  expected. 
Two  types  of  contour  lines  exist  -  “C”  and  “D”.  Note 
that  the  switch  from  the  one  type  to  the  other  occurs 
at  roughly  the  pressure  ratio  value  associated  with  the 
relevant  shock  wave  causing  the  pressure  rise.  The  up¬ 
stream  influence  at  the  line  of  symmetry  is  ~  ISj  for 
Case  F  and  ~  12.Sd  for  Case  G.  For  the  stronger  (8  x  8) 
interaction,  at  the  centerline  near  the  trailing  edge  of 
the  domain,  the  pressure  rise  is  observed  to  reach  a  uni¬ 
form  rate.  As  explained  later,  the  flow  solution  in  this 
region  exhibits  strange  behavior. 

Fig.  30  shows  computed  surface  streamlines  for 
both  cases.  Case  F  exhibits  a  regular  pattern  devoid 
of  any  line  of  coalescence  despite  the  fact  that  the  total 
pressure  rise  after  two  shocks  is  roughly  equal  to  that 
achieved  with  a  single  8  deg  fin.  Outside  of  the  region 
close  to  the  fin  leading  edge,  the  minimum  and  maxi¬ 
mum  flow  turning  angles  are  —3.7  and  17.9  degrees  re¬ 
spectively.  The  surface  streamlines  for  Case  G  form  an 
interesting  pattern.  Near  the  fin  leading  edge,  surface 
turning  angle  roughly  equal  the  shock  angle  resulting 
in  a  dense  region  of  bunched  lines.  Surface  streamlines 
at  and  near  the  line  of  symmetry  do  not  display  signif¬ 
icant  turning.  Upon  approaching  the  line  of  symmetry, 
the  densely  concentrated  lines  straighten  out  and  move 
parallel  to  the  centerline.  Near  the  trailing  edge  of  the 
domain,  the  two  sets  of  lines  &om  either  side  of  the  line 
of  symmetry  form  a  node-like  structure.  Due  to  the  fi¬ 
nite  sise  of  the  grid,  it  is  not  clear  if  there  are  one  or 
two  nodes  (slightly  displaced  from  the  line  of  symme¬ 
try).  We  believe  at  this  time,  that  this  structure  is  a 
numerical  artifact.  The  overall  accuracy  of  the  solution 
near  the  trailing  edge  of  the  domain  is  questionable  (see 
below)  since  the  sublayer  approximation  assumes  uni¬ 
form  mass  flux  within  the  sublayer  region  as  imposed  at 
the  wall  (sero  in  this  case). 

The  computed  fiowfields  are  analysed  with  par¬ 
ticle  traces  as  for  the  single  fin  cases.  The  flowfield 
for  Case  F  exhibits  a  simple  structure  as  summarised 
in  Fig.  31  which  shows  traces  of  particles  released  at 
v/^eo  =  0-1  &nd  y/^ao  =  0.5.  At  the  lower  distance,  the 
particles  near  the  fin  yaw  mote  toward  the  centerline 
than  those  particles  emanating  &om  higher  in  the  up¬ 
stream  boundary  layer  i.e.,  the  spanwise  displacement 


effect  is  larger  near  the  base  of  the  fin.  As  anticipated, 
this  interaction  does  not  exhibit  separation. 

Fig.  32  shows  traces  of  particles  released  for 
Case  G  inside  the  sublayer  at  y/5go  =  0.0001.  The  par¬ 
ticles  follow  the  general  behavior  of  the  surface  stream¬ 
lines  and  form  an  unphysical  stagnation  point  on  the 
plate  at  the  location  corresponding  to  the  node.  Fig.  33 
shows  a  sequence  of  particles  released  at  various  heights 
in  the  boundary  layer.  Particles  released  close  to  the 
plate  but  outside  the  sublayer  (y/^oo  =  0.01)  also  foUow 
the  general  pattern  of  the  surface  streamlines  up  to  the 
node  at  which  point  they  rise  and  move  parallel  to  the 
center  line.  Particles  released  in  the  vicinity  of  the  line 
of  symmetry  at  various  heights  move  parallel  to  it  (as  ex¬ 
pected  &om  the  imposed  boundary  conditions)  but  rise 
up  and  away  from  the  plate  to  form  an  arch-like  trace  at 
the  downstream  end  of  the  domain.  Particles  released 
near  the  fin  at  small  values  of  y/^oo  pitch  toward  the 
plate  and  yaw  toward  the  line  of  symmetry  straighten¬ 
ing  out  as  they  approach  it.  At  higher  distances  from 
the  plate,  particles  released  near  the  fin  pitch  toward 
the  plate  over  the  particles  beneath  yawing  toward  the 
centerline.  The  arch-like  structure  formed  by  the  par¬ 
ticles  released  near  the  centerline  broadens  so  that  the 
overall  fiow  approaches  the  boundary  layer  fiow  (due  to 
the  fin)  expected  at  large  distances  from  the  plate. 

T  Conclusions 

The  effect  of  suction  on  the  supersonic  fiow  past 
the  SFC  is  investigated.  Two  fin  angles  20  deg  and  8  deg 
are  considered,  both  at  Mach  3.  The  former  represents  a 
strong  interaction  for  which  the  existence  of  the  vortical 
structure  has  been  previously  demonstrated.  The  latler 
is  an  incipiently  separated  configuration  for  which  ex¬ 
perimental  data  incorporating  suction  is  available.  For 
the  20  deg  fin,  bleed  is  applied  on  two  regions  on  the 
plate  jointly  spanning  the  region  between  the  line  of  up¬ 
stream  influence  and  the  base  of  the  fin.  For  the  8  deg 
fin,  the  area  of  suction  considered  is  roughly  rectangular 
and  covers  a  strip  on  either  side  of  the  theoretical  shock 
foot  print. 

For  both  interaction  strengths,  sur&ce  turning  an¬ 
gles  are  generally  lower  than  for  the  equivalent  no  bleed 
case.  For  the  stronger  interaction,  the  line  of  coalesence 
persists  for  all  suction  magnitudes  and  areas  considered. 
Since  the  weaker  interaction  does  not  display  a  clear 
line  of  coalescence  in  the  absence  of  bleed,  an  equivalent 
statement  cannot  be  made.  Detailed  examination  of  the 
fiow  indicates  that  the  width  of  the  upstream  compres¬ 
sion  fan  formed  by  the  shock  wave  near  the  plate  is 
significantly  reduced  for  Region  I  bleed,  for  which  up¬ 
stream  influence  is  also  reduced  and  the  rate  of  pressure 
rise  on  the  plate  is  higher  although  the  total  pressure 
rise  remains  roughly  the  same.  Local  skin  friction  val- 
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2.98/7.17 

Table  4:  Theoretical  pressure  rise  alter  shocks  -  Cases  F  and  G 


ues  display  steep  increases  above  no  bleed  values.  For 
the  stronger  interaction,  the  vortical  structure  persists. 
The  effect  of  suction  is  to  ingest  particles  approaching 
the  bleed  area  without  affecting  the  overall  flow  fleld  in 
any  signiflcant  fashion. 

The  flow  past  the  symmetric  double  flu  configu¬ 
ration  is  computed  for  a  weak  and  a  strong  interaction. 
The  weak  interaction  (4  x  4  at  Mach  3)  displays  a  rel¬ 
atively  simple  flow  structure  without  the  existence  of 
separation.  The  flow  past  the  8x8  SDFC  exhibits  sep¬ 
aration  upstream  of  the  primary  shock  with  signiflcant 
flow  motion  away  from  the  plate.  Surface  flow  turning 
angles  exceed  the  inviscid  shock  angle.  Flow  traces  in¬ 
dicate  that  particles  released  at  the  same  height  at  the 
upstream  end  of  the  domain  form  an  arch-like  structure 
at  the  downstream  end  with  the  flow  near  the  plane  of 
symmetry  moving  furthest  away  from  the  plate. 
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•'-FLAT  PLATE  Figure  4:  Symmetric  Double  Fin  Configuration  (SDFC) 
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Figure  1:  Sharp  Fin  Configuration  (SFC) 


Figure  2:  Mean  Flowfield  Structure  for  SFC 
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Figure  3:  Bleed  Regions  -  20"  SFC 


Figure  6:  Experimental  Locations  of  Shapey  and  Bog- 
donoff  -  20°  SFC 
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Figure  W:  Particle  tracea  ia  rerezse  —  Cue  C  (20*  SFC) 
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Figure  18:  Comparison  of  turlkce  ra.w  with 
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Jgure  15:  Particle  traces  -  y/S  =  0.6  -  20 


20  deg  SFC 


Bleed  Region 
Figure  19;  Location  of  experimental  s 


surveys  -  8  deg  SFC 


0.5  o.a 


Legend: 

^ase  D 

. Caafi.E 

O  Sin  1 
&.  =  Stn  Z 
+  =  Stn  3 
X  =  Stn  4 
o  =  Stn  5 
7  =  Stn  6 
B  =  Stn  7 
X  =  Stn  8 
♦  =  Stn  9 


na  Aagi*  ■  •  D«g. 


(b)  Caae  E 

Figure  26:  Particle  tracea  -  yfS  %  0.1  (8*  SFC) 
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Figure  24:  Surface  itreamlincs  -  Caae  E  (8*  SFC) 
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(a)  Caae  D 

Figure  27:  Particle  traces  -  y/S  =  OJ  (8*  SFC) 

(b)  Caae  E 

Figure  25:  Particle  traces  -  y/^  =  0.0001  (8*  SFC) 
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(a)  Case  D 


(a)  Case  D 


(b)  Case  E 

Figure  28:  Particle  traces  -  y/S  =  1.0  (8"  SFC) 
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Figure  31:  Particle  traces  —  Cate  F 
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Figure  32;  Particle  traces  in  sublayer  —  Case  G 
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Figure  33:  Particle  traces  —  Case  G 
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Viscous/Inviscid  Effects  in  3-D  Shock  Wave- 
Turbulent  Boundary  Layer  Interactions 

N.NkrkjaniwBiiu^  D.D.Knight) 

Depurtment  of  Mechanical  and  Aetoapaee  Engineering 
Rutgeri  UniTerrity,  Piacataway,  NJ  088S5. 


AhitnO 

A  detailed  analyeia  of  Tiacona  and  inmeid  eflecta  in 
S-D  ahock  ware  •  tarbnlent  bonndaip  lajer  intcractiona 
haa  been  carried  out,  femaaing  on  tbe  (aiA)=  (24*, 40*) 
awept  compceaaioa  comer  eonSgiiratian  at  itea.  =8.1x10* 
and  tbe  «  =20*  abarp  fin  configuration  at  Jlea.  =8.0x10*. 
Two  fiowfidda,  computed  naing  diflerent  tnrbnleace  mod- 
cia,  are  analjraed  for  each  eonfigoration.  Tbe  mean  H- 
netie  cnergp  equation  ia  integrated  along  atreamliaea  orig¬ 
inating  at  varioua  j/iia  loeationa  within  and  outaide  tbe 
boundary  layer  and  a  quantitatirc  eompariaon  ia  made  of 
tbe  contribution  to  tbe  mean  kinetic  energy  by  inriacid 
and  eiaeona  terma  ia  the  equation.  Alao,  tbe  terma  are 
quantitatively  compared  at  pointa  on  y  •  a  planea  (normal 
to  tbe  atreamwiae  directioa)  Cor  the  abarp  fin  interaetiem 
and  at  poanta  lying  on  f  •  a  planea  (along  tbe  atreamwiae 
direction)  Cor  the  ewept  compreaaion  comer  interaction. 
Tbe  reaulta  indicate  that  tbe  iateractiona  Cor  both  geome- 
triea  are  governed  primarily  by  inviacid  (preaaure)  efTeeta. 
Viacoua  infiueace,  although  not  of  tbe  aame  magnitude  aa 
the  inviacid  eflecta,  ia  found  to  be  aignificant  in  regioua 
botb  near  and  away  from  tbe  aur&ce  of  tbe  geometry  and 
not  reatricted  to  tbe  *viaooua  aublayer*  re|pon  within  tbe 
boundary  layer. 

Introdnction 

Tbe  problem  of  S-D  abock  wave  -  turbulent  bound¬ 
ary  layer  iateractiona  (  S-D  turbulent  iateractiona  )  oe- 
cura  ia  a  variety  of  bigb  Mach  number  fiowa.  Some  of 
tbe  practical  applicatioaa  in  which  tbu  phenomenon  ia 
encountered  include  high  apeed  aircraft  inleta  and  com- 
preaaon,  wing-body  jnnetnrca  and  control  aurCacea.  The 
eompleiity  of  aucb  viaeona-iaviaad  iateractiona  ia  doc  to 
non-liaearity,  eomprearibility  eflecta  and  high  gradieata  of 
flow  variaUea,  bcaidea  tbe  difficulty  ariaing  from  three- 
dimenaionality  and  turbulence.  The  iateractiona  arc  gov¬ 
erned  ia  each  eaae  by  tbe  freeatream  Mach  number  iimt 
tbe  Beynotda  number  Jfea.  baaed  upon  tbe  boundary  layer 
tbickneca  the  wall  thermal  eonditioua  and  the  geome¬ 
try  generating  tbe  interaction. 
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Significant  progreaa  baa  been  made  ia  tbe  underatandiag 
of  S-D  turbulent  intcractiona  ia  tbe  pact  Cew  yeara.  Inter- 
actiona  generated  by  dimenaionleea  geometrica  like  abarp 
fina  (Fig.l)  and  awept  compreaaion  comen  (Fig.2)  have 
been  atudied  extenaivriy  naing  experimental  teebniquea, 
numerical  flow  aimulation  and  analytical  metboda. 

Recent  experimental  inveatigationa  of  S-D  turbulent 
iateractiona  for  rimplified  geometrica  like  abarp  and  Uunt 
fina  and  awept  compreaaion  eomen  have  been  carried  out 
by  a  number  of  rcaeareben.  Experimenta  with 

aurfoce  and  boundary  layer  meaanremmta  have  bem  con¬ 
ducted  by  McCabe  (1),  Law  [2],  Xorkegi  [S],  Peake  (4), 
Kuaaop  ct  a/  [S|,  Kubota  and  StoDery  [8],  Otkarn  ct  W 
(7,8|,  Dolling  and  Bogdonoff  [8,10),  McClure  and  Dolling 
(11)  and  SetUca  ct  al  [12,18,14,15,18]. 

Analytical  atudiea  of  S-D  intcractiona  have  been  made 
by  Myring  |17|,  Stalker  [18,18),  Inger  [20,21,22]  and  by 
Mee  and  Stalker  [23]. 

Numerical  aimulationa  uaing  tbe  3-D  mean  com- 
preaaible  Reynolda  averaged  Navicr  -  Stokea  cquationa  have 
focuaaed  on  aimilar  dimenaioaleaa  geometrica.  Swept  com¬ 
preaaion  comer  and  abarp  fin  generated  iateractiona  have 
been  atudied  by  Hung  and  MacCormack  [24,25],  Horatmnn 
and  Hung  [27],  Horatman  [28,28]  and  Knight  [30,31,32,33]. 
The  computationa  of  Knight  employ  tbe  Baldwin-Lomax 
algebraic  turbulent  eddy  viacoaity  mcxlel  [30]  and  Iboae 
of  Horatman  incorporate  tbe  Jonea-Launder  (k  -  c)  two 
equation  model  [37]  or  Cebici-Smith  model  [38].  Theae 
computationa  have  been  validated  by  eompariaon  with  ex- 
perimenta  and  have  been  found  to  be  ia  good  agreement 
with  meaaured  vuluea  of  pitot  and  atatie  preaaure  profilea 
ia  tbe  boundary  layer,  aurfoce  preaaure  and  yaw  aaglea 
[34,35]. 

Tbe  iateractiona  are  found  to  be  dominated  by  a  large 
vortical  atructurc  (Figa.  5  and  6),  ia  agreement  with  t^ 
flowfield  modcla  propoacd  by  Token  [40]  and  and  Kubota 
and  StoUery  ]6].  Bcaidea,  tbe  computationa  of  Knight  and 
Horatman,  although  employing  two  diflerent  turbulence 
modcla,  abow  cloae  agreement  for  botb  abarp  fin  and  awept 
compreaaion  comer  iateractiona,  diaplaying  modeat  difler- 
eacea  only  ia  a  amall  fraction  of  tbe  boundary  layer  ad¬ 
jacent  to  the  aurfoce.  The  ifuenaitivily  »f  tke  flemfitU  to 
Ike  iurkuience  models  employed  m  tke  compulatiofu  auf- 
ftsU  lAal  tke  flowfield  is  primarily  rotational  and  inviacid 
witk  viacoua  influence  keiny  limited  to  rcyioru  very  close 
to  tke  wall  or  surface. 
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Ib  order  to  verify  this  conjectiire  ead  to  further  under- 
ttud  the  dynamica  of  the  isteractioot,  the  present  study 
aims  at  making  a  fuantitaiise  comparison  of  invisdd  and 
viscous  influence  in  the  flow.  Attention  has  been  focussed 
on  the  terms  in  the  mean  kinetic  energy  equation  (chosen 
tor  simplicity  and  aa  a  first  step  in  this  quantitative  anal¬ 
ysis)  and  a  computer  code  has  been  developed  to  analyse 
the  contributions  to  the  rate  of  change  of  mean  kinetic 
energy  hy  both  invisdd  and  viscous  terms,  along  stream¬ 
lines  in  the  flow.  In  particular,  this  code  has  been  used 
to  analyse  the  3-0  turbulent  interactions  generated  by  a 
swept  compression  comer  and  a  qualitatively  similar  in¬ 
teraction  generated  by  a  sharp  fln  geometry.  Both  geome¬ 
tries  chosen  are  'dimensionlcas'  and  the  interactions  in  the 
two  eases  involve  similar  incoming  flow  characteristics  and 
pressttie  rise  (shock  strengths). 

Computation  of  Flowfleldi 

The  flowfidds  for  the  sharp  fln  and  swept  eompres- 
sion  comer  interactions  analysed  in  the  present  study  have 
been  computed  by  Knight  [30,31]  and  Gaitondc  [41]  and 
separatdy  by  Rorstman  [28,20].  The  mean  compress¬ 
ible  Navicr-  Stokes  equations  coupled  with  the  Baldwin- 
Lomax  algebraic  turbulent  eddy  viscosity  modd  (Knight 
and  Gaitonde)  or  the  Jones- Launder  (k  —  r)  two  equation 
modd  (Rorstman),  form  the  basis  of  these  computations. 
The  code  employed  by  Knight  uses  an  efficient  hybrid 
cxplidt-implieit  algorithm  [30]  for  sdution  of  the  govern¬ 
ing  equations.  The  implicit  algorithm  is  the  Box  Scheme 
of  KeUer  [30],  which  is  applied  to  the  asymptotic  form 
of  the  Navier-Stokes  equations  in  the  very  thin  ‘viscous- 
sublayer’  region  of  the  boundary  layer.  The  flowfidd  in 
the  region  outside  the  sublayer  is  computed  using  the  ex¬ 
plicit  finite-dificrenee  algorithm  of  kfaeCormack  [25,26]. 
The  computations  incorporate  a  body-fitted  coordinate 
transfocmation  to  permit  the  treatment  of  flow  regions  of 
arbitrary  shapes  and  the  code  has  been  vectorised  for  d- 
fldent  computation  on  vector  processing  sopereompuiers 
like  the  CYBEIl-208  and  ETA- 10.  The  code  employed  by 
Rorstman  uses  an  explicit  finite-diflerence  scheme  for  the 
computations. 

The  present  study  involves  the  examination  of  the  Col- 
lowing  flowfidds  in  detail: 

1.  (arl)a(24*,40')  Swept  comer  at  Rei.  =8.1x10*  , 
Jfia=3  and  fa>  =  0.5  inch. 

(Flowfield  No.l:  Computed  using  the  Bddwin-Lomax 
modd,  and  Flowfidd  No.2:  Computed  "=«g  the  Jones- 
Launder  modd). 

2.  a  =20"  Sharp  fin  at  Res.  =0.0x10*,  Af,.=3  and  <.  = 
0.5  inch. 

(Flowfidd  No.l:  Computed  using  the  Bddwin-Lomax 
modd,  and  Flowfidd  No.2:  Computed  using  the  Jones- 
Launder  modd). 


Governing  Equation  for  Mean 
Kinetic  Energy 

The  mean  kinetic  energy  equation  dong  a  streamline 
may  be  written  aa  fallows: 


Rere,  |ujUi  is  the  mean  kinetic  energy  per  unit  mass. 
The  D/Dt  operator  denotes  the  Particle  derivative  or  the 
Substantid  Derivative,  p  denotes  the  density  of  the  fluid, 
p  U  the  pressure  and  (i=l,2,3)  denotes  the  coordinates 
X,  y  and  s. 

r^  is  the  stress  tensor  (laminar  +  turbulent)  and  is  given 
by  the  rdation: 


m  - -fo*  +  «)V  •  fy  +  (p  +  .)(^ -i- ^) 

where, 

tii=  1  for  i=j  ,  and  =  0  for  i  j 

l.j=lA2 


c  is  the  algebrde  turbulent  eddy  viscosity  given  by  the 
turbulence  modd  employed,  is  is  the  velocity  vector  and  p 
is  the  molecular  dynamic  viscosity  obtained  Suther¬ 
land's  law. 

The  equations  above  are  written  in  the  summa¬ 

tion  convention. 

lolcgration  of  the  mean  kinetic  energy  equation  yidds: 


-MiUi  =  ju,m  + 


/  *<  ^  **  I 

J  ~ pBziV  J  pflx,  V 


where, 

{waW,  denotes  the  mean  kinetic  energy  at  the  starting 
point  of  streamline  for  the  fluid  particle.  'V'  denotes  the 
magnitude  of  the  resultant  vdodty  vector  at  the  location 
of  the  fluid  partide  and  ‘ds’  is  the  ii«6m»Tfiiml  ardength 
dong  the  streamline. 

The  rate  of  change  of  mean  kinetic  energy  per  unit  =«»■« 
is  aflected  by  two  distinct  terms: 

The  first  term  — (uj/ p)dp/dzi  denotes  the  invisdd  or  pres¬ 
sure  efiect  and  the  second  term  (u</p)flr^/aa^  denotes  the 
combined  viscous  and  turbulent  stress  efleet.  By  following 
a  partide  dong  a  streamline  and  simultaneoudy  integrat¬ 
ing  the  mean  kiaetie  energy  equation  along  the  path,  the 
magnitudes  of  the  invisdd  contribution  and  viscous  con¬ 
tribution  to  kiaetie  energy  are  qnaatiutivdy  compared. 


It  is  also  useful  to  eondder  the  following  two  groups  of 
terms: 


1. 

2. 


I 


I 


Bw  fln  ds 


The  first  group  denotes  the  mean  h'nelic  cneryy  ok- 
tained  ky  inefudinf  kotk  ike  mtucid(prt$$urt)  and  the  vii- 
cous  inieyrsndi  and  the  second  group  denotes  the  mean 
kinciie  eneryy  okisined  ky  tnc/udiny  only  ike  pressure  inie- 
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ffnd.  Th*  groups  1  ud  2  we  lunultMeousIy  computed 
•long  the  pnth  of  •  pwticle  (along  the  •treamline)  and 
competed  with  each  other.  A  good  agreement  between  the 
two  result*  indicates  that  the  dTnamics  along  the  streanr- 
Une  are  goeemed  prindpallp  by  the  ineueid(pretsure)  ef¬ 
fects,  while  a  ugnificaat  disagreement  between  the  two 
indicates  that  &olA  prsMure  and  tucous  effteU  influence 
the  flow  along  the  streamline.  Various  streamlines  (far 
particles  starting  at  diSerent  points  in  the  flow)  ate  in- 
reatigated  in  this  manner  and  the  results  analysed  to  ver¬ 
ify  the  original  conjecture  about  the  prinapally  rotational 
and  ineiadd  nature  of  the  flowfldd. 

Computational  Khenw  for 
tha  machanical  cnargy  equation 

Both  the  swept  compression  comer  and  the  sharp  fin 
geometries  ineolve  aon-nniCatm  grids  in  the  ‘phydcal*  do¬ 
mains  of  the  flow.  As  mentioned  earlier,  in  computing  the 
flowflelda,  a  body-fitted  coordinate  transformation  is  used 
to  permit  treatment  of  these  non-uniform  regions.  In  a 
manner,  the  eraluation  of  the  path  integral  (of  the 
mean  kinetic  energy  equation)  involTes  the  transformation 
of  the  aon-uniibrm  ‘physical’  grid,  denoted  by  the  (x,y,s) 
coordinate  system  to  a  uniform  'computational*  grid  de¬ 
noted  by  the  ((,i|,()  coordinate  system.  The  derivatire 
terms  in  the  equation  are  then  ^proamated  by  second- 
order  accurate  flnite-diflereaces.  The  partide  trsdng  is 
accomplished  by  an  algorithm  developed  by  Raufia'(42}, 
which  involves  location  of  the  partide  within  a  partidar 
‘cell’  of  the  grid,  interpolation  of  the  vdodty  terms  at 
the  particle  position  and  ‘movement’  of  the  partide  using 
the  interpolated  vdodties.  This  partide  tracing  code  has 
been  used  to  study  the  flowfidd  structure  of  the  (ci,A)s 
(24*,40*)  swept  compression  comer  interaction  (3Sj.  The 
mean  kinetic  energy  equation  is  integrated  simultaneaudy 
with  each  movement  of  the  particle.  The  iavisdd(pressure) 
and  viscous  (laminar-t-turbulcat}  contributiona  to  the  rate 
of  change  of  mean  kinetie  energy  are  computed  at  the  pw- 
tide  podtion  in  the  following  steps: 

1.  Location  of  the  partide  in  the  uniform  (computational) 
grid  (using  the  partide  tracing  algorithm) 

2.  Evaluation  of  the  invisdd  and  viscous  terms  of  the 
mean  kinetie  energy  equation  at  the  ‘cell’  comers  using 
flnite-diflSrrences. 

S.  Linear  interpolation  of  these  vdues  at  the  partide  po¬ 
sition  and 

4.  Integration  of  the  invisdd  and  viscous  terms  (at  the 
partide  podtion)  udng  simple  trapesoidd  method. 


(a,A)=:(24*,40*)  Swept  Compression  Corner 
/fe«.  =  8.1x10* 

The  swept  comptesdon  comer  geometry  and  the  com- 
putationd  domain  are  shown  in  Pip.2  and  4.  The  shock 
generated  by  the  presrace  of  the  swept  comer  interacts 
with  the  incoming  turbulent  boundary  layer  on  the  flat 
plate  and  causes  sepwation  of  the  flow.  The  computation 
of  the  flowfldd  incorporates  the  same  incoming  flow  condi¬ 
tions  as  the  experiments  conducted  far  this  configuration 
at  the  high  Reyndda  number  wind  tnnnd  at  the  Prince¬ 
ton  Gas  Oyaamics  Laboratory.  The  experiments  involved 
a  Mach  number  of  3,  a  totd  pressure  of  690  kPa  ±  1.3%, 
a  totd  temperature  of  266.4  deg  K  ±  1.6%  ,  near  adia¬ 
batic  surface  temperature  and  an  incoming  2-0  turbulent 
boundary  layer  which  has  been  studied  and  observed  to 
satisfy  the  Law  at  the  Wall  and  Wake  [43|.  The  molecular 
and  turbulent  Prandtl  numbers  are  0.73  and  0.9,  tespee- 
tivdy.  The  incoming  tnrbulmt  bounday  layer  is  a  ‘thick’ 
layer  of  a  height  s  0.3  inch.  Here,  refers  to  the 
height  of  the  boundary  layer  at  the  comer  line  (where  the 
compression  ramp  meets  the  flat  plate),  in  the  absence  of 
the  ramp. 

Flowfldd  Wo.l:  fle/dwin-Iomaa  Ihrfu/ewce  Hotel 

In  order  to  study  individud  partide*  in  the  flow,  stream¬ 
lines  were  traced  by  rdeaaing  partides  at  fixed  j/Sm  loca¬ 
tions  near  the  upstream  end  of  the  flowfldd.  10  partide* 
were  tdeaaed  at  every  y/f.  location,  spaced  evenly  in  the 
s-direetion,  between  the  minimum  and  maximum  s  vdues 
of  the  grid  at  that  location.  Partide  traces  for  ys0.02f,„ 
0.0Sf«,  0.2tm  sad  0.Sf«  are  shosrn  in  Pig.7.  The  parii- 
de*  rdeased  at  y=0.02f«  and  y=0.056a  sre  observed  to 
get  entrained  into  the  vortied  straeture  at  the  base  of  the 
compressios  comer  and  those  starting  at  ys  fi.Hm  and 
y=0.Sf«  continue  over  the  vortex  and  up  the  compres¬ 
sion  ramp. 

The  streamlines  may  be  numbered  1  through  10,  start¬ 
ing  from  the  one  nearest  to  the  symmetry  boundary  of  the 
domain.  The  /(/U  streamline  among  these  is  chosen  as  a 
‘typied’  streamline  tor  study,  for  each  of  the  four  y/f« 
(starting)  locations.  In  Figs.  8-11,  the  kinetic  energy 
evduated  by  interpolating  the  vdodties  ficom  the  com¬ 
puted  flowfldd  and  that  obtained  by  computing  the  path 
integrd  (right  hand  side  of  the  integrated  mean  kinetie 
energy  equation)  are  shown  for  each  y/f«  (starting)  case. 
The  results  for  the  ys  0.026^  ,y=0.056a  and  y=0.66« 
cases  have  been  presented  in  an  earlier  paper  [38],  but 
are  described  again  here  for  comparison  with  the  results 
for  the  flowfldd  obtdned  using  the  different  turbulence 
modd  (Flowfldd  No.2).  The  kinetic  energy  profiles  are 
in  dose  agreement,  and  the  differences  in  the  two  results 
is  attributable  to  the  truncation  errors  involved  in  the 
schemes  employed  in  computing  both  the  path  integrd  as 
well  as  the  flowfldd  itself.  For  the  streamhnes  originating 
at  y=0.26M  and  y=0.3f«,  however,  the  difference  in  the 
two  vdues  of  kinetie  energy  is  spproximatdy  16%  in  each 


3 


eaac.  The  reaion  tar  thia  error  ia  being  ineeatigated.  Tbe 
kinetic  energy,  in  eaeb  of  the  caaca,  remaina  unaffected 
till  the  itart  of  the  interaction,  when  the  adeerae  preaaure 
gradient  encountered  in  the  region  of  the  ihoch  coupled 
anth  eiacoua  influence  cauaea  a  aharp  decline  in  the  energy 
of  the  particle.  The  drop  in  energy  itopa  at  the  end  of  the 
interaction,  after  which  it  rcachea  a  neat  aaymptotie  ealue. 

For  the  came  atreamUnea  diacuaaed  aboac,  the  mean 
kinetic  energy  obtained  by  including  both  the  preaaure 
and  the  eiaeoua  integranda  ia  compared  to  that  obtained 
by  the  path  integral  iaaoleing  only  the  prcature  integrand, 
(i.e)  the  term 

i-i-t  I—  +  +/ 

ia  compared  to  the  term  |uin< 

The  flrat  term  ia  referred  to  aa  K.E  -Ptik  IiUtgrat,  in 
the  plota  and  the  aecond  term,  aa  E.t  -Path  InUftmt,  /la- 
fuccd  /wieyreiad  For  the  ys0.02la  caae,  the  two  ealnaa 
differ  eonaiderably,  indicating  that  kett  maiactd  and  eia* 
c««a  ejecta  are  ri^ijSeanl  there  rtreamimer.  Similar 
reanita  ate  obeerred  fee  the  other  nine  atreamJinea  orig* 
mating  at  y=0.02dm.  For  the  ysO.OSds  eaae,  the  two 
raluea  are  reaaonably  eloee  and  it  can  be  inferred  that  the 
prearere  (inwUcid)  in/lucnce  ir  deminant  Car  the  atream- 
line.  Similar  reaulta  are  again  obaerred  Cor  the  other 
ctreamlinea  originating  from  the  came  y-  loea_on.  It  ia  to 
be  noted  here,  that  in  both  the  y=0.02fa  and  y=O.OSds 
caaea,  the  particlea  get  trapped  deep  within  the  Torticai 
atructure.  Yet,  riacoua  influence  ia  more  aigniflcaat  Cor  the 
former  than  the  latter  caae.  The  plota  Cor  the  ysO.Sda  aad 
ysO.Sda  ctreamlinea  ahow  an  inoreaaiim  dominance  of  the 
ineiadd  effect  on  the  mean  kinetic  energy.  Aa  remarked 
beCore,  theae  particles  do  not  get  trapped  within  the  rot- 
ticai  atructure  of  the  interaction,  hut  continue  over  it  and 
along  the  ramp.  The  remaining  nine  ctreamlinea  ia  each 
eaae  exhibit  aimilar  behavior  aa  the  ‘typical’  atreamliae 
examined  Cor  the  two  caaea. 

The  relatiTe  importance  of  the  viscous  effect  Cor 
the  streamlines  originating  at  y^O.OTt^  and  ysO.OSd«  ia 
also  evident  from  the  proflle  of  the  kanetic  energy  obtained 
using  only  the  viscous  integrand  (reCeied  to  aa  JT.ff  -FotA 
Integral,  Vucour  Integrand)  ia  the  plots. 


A  ajmila*  analysU  is  carried  out  Cor  thia  flowfield  and 
the  mean  kinetic  energy  evaluated  Cor  streamlines  orig¬ 
inating  at  frsO.OTS^,  0.0M«,  0,2Sm  and  Again, 

the  flfth  streamline  from  the  symmetry  boundary  of  the 
domain  ia  chosen  as  the  ‘typical’  atreamliae  Cor  study. 
Figs.13  and  Ifl  show  the  ctreamlinea  (or  the  particles 
originating  at  the  Cour  different  y  -locations.  The  atream- 
lines  are  similar  to  that  obtained  Cor  the  different  cases  Cor 
Flowfleld  No.l.  Figs.lS-14  and  ld-17  ahow  a  comparison 
between  the  kinetic  energy  obtained  by  a)  interpolating 
the  flowfleld,  b)  by  the  path  integral,  c)  by  the  path  in¬ 
tegral  obtained  with  only  the  inviacid  integrand  and  d) 
by  the  path  integral  obtained  with  only  the  viscous  inte¬ 


grand.  The  kinetic  energy  values  obtained  by  flowfield  in¬ 
terpolation  and  by  tbe  total  path  integral,  although  close, 
differ  by  about  10-157(  on  a-^age.  This  error  may  be 
attributed  to  the  coarseness  of  the  grid  and  the  resulting 
truncation  errors  in  the  computation  of  both  the  flowfleld 
and  the  path  integral.  Consequently,  no  definite  conclu¬ 
sion  can  be  reached  about  the  nature  of  the  flow  from 
these  reaulta.  It  can  however  be  noted  that  the  path  inte¬ 
gral  obtained  using  the  inviacid  term  alone  approximates 
reaaonably  well,  the  total  path  integral  in  each  of  the  Caur 
cases.  Also,  the  viscous  influence  ia  particulariy  aignifl- 
cant  Cor  the  atreamliae  originating  at  y  =  0.2Sm- 

30*  Sharp  flat  flea.  =  ^0x10* 

The  aharp  fin  geometry  and  the  computational  do¬ 
main,  arc  shown  ia  Figs.  1  and  S.  Bcre  again,  the  gener¬ 
ated  shock  wave  intcracU  with  a  ‘thick’  ineoaniag  bound¬ 
ary  layer  on  the  flat  plate  (<•»  0.S  inch),  where  f.  is  the 
height  of  the  boundary  layer  at  the  location  of  the  apex 
of  the  fin,  ia  the  absence  of  the  fin. 

The  computations  of  the  flowfield  incorporate  the  in¬ 
flow  conditions  Cor  the  experiments  on  the  20^  aharp  fin 
caae  at  the  Princeton  Gas  Dynamiea  Laboratory.  The 
experiments  involve  nominal  freeatream  Mach  number  of 
3.03,  a  total  pressure  of  830  kPa  ^  1%  ,  a  total  temper¬ 
ature  of  251  deg  K  ±  i%,  near  adiabatic  wall  conditions 
and  a  2-0  ‘thick’  incoming  equilibrium  boundary  layer 
which  has  been  thoroughly  investigated  and  observed  to 
fit  the  Law  of  WaU  and  Wake. 


An  analysis  to  that  for  the  swept  compression 

comer  caae  is  carried  out  Cor  thia  interaction.  10  parti¬ 
cles  spaced  evenly  between  the  minimum  and  maximum 
's’  values  of  tbe  domain  upstream  of  the  leading  edge  of 
the  fin  are  released  at  different  y/f«  locations  near  the 
upstream  boundary.  Again,  the  'fiftk’  streamline  is  cho¬ 
sen  (or  study  (or  each  of  the  y-levela.  The  streamlines  fer 
particles  released  at  y=0.02fv  ,y=0.0Sf,s,  y=0.2fM  and 
are  shown  in  Flg.18.  In  all  the  (bur  caaea,  the 
streamlines  are  observed  to  get  wrapped  up  into  the  vor¬ 
tical  structure  emanating  from  near  the  leading  edge  of 
the  fin.  Figa.lB-22  show  the  comparison  between  a)  the 
mean  kinetie  energy  obtained  by  interpolati«»  of  the  com¬ 
puted  flowfield  b)  the  path  integral  (or  the  streamline  e) 
the  path  integral  involving  only  ine  inviacid  term  and  d) 
tbe  path  integral  involving  only  the  viscous  term,  (or  each 
eaae.  The  interpolated  and  integrated  values  of  kinetie  en¬ 
ergy  agree  quite  closely  in  the  (bur  eases,  and  differences 
in  the  two  values  are  attributable  to  the  truncation  errors 
in  the  numerical  schemes.  The  y=0.02fM  ease  shows  a 
drop  in  the  kinetie  energy  at  the  start  of  the  interaction 
but  a  subsequent  rise  ia  the  energy  to  a  level  above  the 
value  at  the  upstream  (starting)  position  of  the  particle, 
i.e.,  the  particle  ia  ‘accelerated’  by  the  influence  of  (he  in¬ 
viacid/ viscous  terms.  The  y=0.05fa  exhibits  a  similar  be¬ 
havior  of  acceleration  of  the  particle  along  the  streamline. 
For  the  g=0.2S„  and  j=0.5i^  cases,  the  mean  kinetie  en¬ 
ergy  drops  at  the  start  of  the  interaction  and  reaches  a 
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near  nyimptotie  value  thereafter. 

It  is  also  observed  that  the  results  obtained  by  comput¬ 
ing  the  Kill  path  integral  (b)  and  the  path  integral  involv¬ 
ing  only  the  invisdd  integrand  (e)  difer  ngnificantly  tot 
the  streamlines  starting  at  y  s0.02i_  and  ysO.OS4«.  Also, 
similar  behavior  is  exhibited  by  the  other  nine  stream¬ 
lines  in  the  two  eases.  It  can  therefore  be  concluded  that 
both  insweid  and  suceus  effteU  art  sifnifieant  tat  these 
streamlines. 

For  the  y=0.2da  and  y=0.Sd«  cases,  the  path  integral 
involving  the  pressure  term  alone,  is  a  reasonable  approx¬ 
imation  to  the  mean  kinetie  energy  of  the  particle,  indi¬ 
cating  that  lAe  /tarn  eisny  ttese  streamlines  it  fosemed 
prineipsll^  hp  invuetd  ffTtttan)  tfftcU.  Similar  results 
are  obtained  tat  the  other  nine  streamlines  originating  at 
these  two  jjtm  locations. 

Figi.23  and  24  show  the  magnitndea  of  the  invisdd 
(pressure)  contribution  to  mean  kinetie  energy  (i.e), 

/— (•s/p)ftp/fts<(|i)  and  the  viscous  contribution  to  mean 
kinetic  energy  (i.e),  /(n</p)^s/ftai(lF)  ^  ‘fifth*  stream¬ 
lines  starting  at  the  levels  y=0.02fa  and  y=O.OSfa-  The 
hdght  of  the  particles  (from  the  surface  of  the  fiat  plate) 
are  also  shown  as  a  function  of  arclength.  The  invisdd 
and  viscous  terms  have  been  normalised  by  the  value  of 
the  mean  ktnetie  energy  at  the  upstream  (starting)  loca¬ 
tions  of  the  partides.  It  is  observed  that  viscous  influence 
continues  to  remain  significant  even  as  the  partide  moves 
sway  from  the  wall,  contrary  to  the  original  hypotheds 
that  a  ‘decrease’  in  the  viacous  efCect  should  be  observed 
as  the  distance  of  the  partide  from  the  wail  increases. 

Tkit  MieaUt  tkai  Iktrt  esisl  reyseiu  in  the  flam,  apart 
from  rtfiant  adjactnt  ta  the  piste  suf/bce,  where  the  in- 
tersetien  it  tifnificanUp  titeoat  in  nsturc. 

Flowfldd  t(o.2;  Jonet-Laanier  UodtL 

Fig.2fi  shows  the  streamlines  tot  partides  rdeased  at 
ys  0.02f_,  O.Ofida,  0.2fa  and  0.9d«.  The  streamlines 
are  similar  to  that  obtained  tot  the  corresponding  cases 
tot  Flowfldd  No.l.  As  befiwe,  the  ‘fifth*  streamline  from 
the  symmetry  boondaiy  is  chosen  toe  study  in  each  case. 
Fig.2fi-27  show  the  kinetie  energy  computed  by:  a)  in¬ 
terpolating  the  velodties  from  the  known  flowfldd  values, 

2)  the  path  integral,  3)  the  path  integral  involving  only 
the  invisdd  integrand  and  d)  the  path  integral  invdving 
only  the  viscous  integrand,  (or  the  fiflk  streamline  origi¬ 
nating  at  y=0.2f«  and  0.33,..  The  kinetie  energy  values 
obtained  by  interpolatioa  and  by  the  path  integral  are 
reasonable  doee.  The  kinetie  energy  profile  obtained  by 
employing  the  the  invisdd  term  alone  provides  a  good  ap¬ 
proximation  to  the  value  obtained  by  the  full  path  integral 
in  both  eases.  The  other  nine  streamlines  in  each  case  ex¬ 
hibit  dmilar  behavior  and  it  can  thus  be  eonduded  that 
the  dynamics  of  the  interaction  along  these  streamlines 
are  governed  primarily  by  inaiteiJfprittari)  ejects.  The 
computations  tot  the  y=  0.023.  and  0.053.  eases  did  not 
yidd  a  good  agreement  between  the  interpolated  and  path 
integral  values  owing  to  large  truncation  errors  in  the  nu¬ 
merical  schemes  and  the  results  for  these  streamlines  are 
hence  not  presented  here. 


PlanewUe  examingtion  of  inviicid  and 
viacous  influance  in  the  flow 

Cass  li  30*  Sharp  Fin  Intcrnctioa 

la  order  to  identify  regions  of  invisdd  and  viscous  in¬ 
fluence  in  the  flow,  the  pressure  term  ~{ai/p)dp/dzi  and 
the  viacous  term  (ai/p)9r,i/8aj  are  computed  at  mesh 
points  oa  p  -  t  planes  at  diflereat  ‘x/3.*  locations  in  the 
flow. 

The  following  analysis  refers  to  the  flowfldd  generated 
using  the  Baldwia-Lomax  algebraic  turbulent  eddy  viscos¬ 
ity  modd  (Flowfldd  No.l).  Fig.2fl  shows  the  magnitude  of 
the  pressure  and  viscous  terms  in  the  form  of  contour  plots 
on  the  p  -  s  plane  at  the  location  xs  153.,  dowastresun 
of  the  leading  edge  of  the  fltu  The  minimnm  s/3,  value  of 
0.0  corresponds  to  a  point  0.013.  away  from  the  surface 
of  the  fin  in  the  s-direction.  The  minimum  y/3.  value  of 
0.0  also  corresponds  to  a  point  0.013.  above  the  surface  of 
the  plate.  The  value  of  y  is  2.03..  The  p  and  s 

values  arc  denoted  as  multiples  of  3..  Both  pressure  and 
viscous  term  values,  ~(ai/p)9p/dat  and  (nifp)6tij/daf 
are  normaliaed  by  the  majdmnm(absolnte)  value  of  in- 
viadd/viscous  contribution  terms  on  the  particular  p  —  s 
plane.  Also  as  observed  from  the  plot  of  ‘prcsanie  con¬ 
tours*  at  the  same  location,  the  pressure  rise  due  to  the 
shock  is  distributed  within  the  boundary  layer  and  forms 
a  *A*  shaped  region. 

A  comparison  of  the  contour  plots  of  the  inviscid  and 
viscous  contributions  (A  and  B)  and  pressure  (C)  shows 
that  the  inviscid  influence  is  concentrated  only  towards 
the  freestream  side  of  the  ‘A’  shaped  shock  structure,  al¬ 
though  the  pressure  rise  continues  upto  the  fin  surface. 
The  reyion  af  viscous  inftatnea  is  tifnificantlt  lar/a  as  is 
evident  from  plot  B  and  lies  to  the  left  of  the  region  of 
muiinum  inviscid  influence.  The  reason  tot  the  concen¬ 
tration  of  inviscid  influence  towards  the  freestream  side 
of  the  ‘shock*  region,  becomes  evident  when  the  inviscid 
(pressure)  term  -(a,'/p)dp/dx<  is  examined  more  closely. 
Since  the  term  involves  both  the  directions  and  magni¬ 
tudes  of  velocity  and  the  pressure  gradient,  the  two  may 
be  combining  in  a  way  so  as  to  result  in  an  almost  ‘rero* 
inviscid  effect  in  the  area  between  the  observed  region  of 
concentrated  inviscid  influence  and  the  flu  surface.  Plot  O 
shows  the  direction  of  the  velocity  vectors  on  the  plane  as 
observed  from  the  apex  of  the  fin.  The  resultant  velocities 
at  the  points  on  this  plane  are  obtained  from  the  velocity 
components  tangential  to  a  ray  drawn  from  the  apex  to 
each  grid  point  on  the  plane. 

The  orientation  of  the  velocity  vectors  clearly  indicates 
the  presence  of  the  vortical  structure  and  the  rotational 
nature  of  the  flow.  In  the  region  between  a  ficticious  'cen¬ 
ter*  of  the  vortical  structure  and  the  fin  surface,  the  flow  is 
observed  to  be  going  ‘down*,  towards  the  plate.  It  is  to  be 
noted  that  at  each  mesh  point  on  the  y-s  plane,  a  compo¬ 
nent  of  velocity  along  the  ray  drawn  from  the  apex  of  the 
Bn  to  the  grid  point,  exists,  but  is  not  represented  in  the 
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diagrua.  It  is  conjectuicd  that  the  tens  -(«i/p]dp/d«4 
becomca  nearly  lero  in  thia  region  due  the  (act  that  the 
dircctiona  of  the  ecloaty  eoaponeata  are  normal  to  that 
of  the  corretponding  preaaure  gradient!  in  thia  region. 
Similar  reaulla  are  obtained  (or  the  other  y- a  planea  anal* 
yacd  at  other  a/f^  location!  in  the  flow.  The  region  of 
viacoua  influence  ia  aigniflcaat  in  each  of  the  other  caaea, 
aa  obaerred  earlier  (or  the  z=  15f_  plane  and  alwaya  liea 
to  the  immediate  left  of  the  region  of  concentrated  inria- 
dd  influence  ia  the  flow. 

The  maximum  abaolnte  ealue  of  the  riacoua  term  ia 
appmamately  409(  of  the  abaolute  ealue  of  the 

iaeiadd  effect  (or  poiata  on  the  zw  plane.  Beaidea, 
the  region  of  eiacoua  influence  extenda  upto  a  hdght  of 
2.04^  aboea  the  aurfaee  of  the  plate. 

The  aaalyaia  (or  Flowfleld  No.2  waa  carried  out  ia  a 
manner  aimilar  to  that  Cor  Flowfleld  Ilo.l.  (reaulta  not 
ahowa  here).  It  waa  obaerred  that  although  the  location  of 
the  region  of  iaeiadd  influence  waa  found  to  agree  cloeely 
with  that  Car  Flowfleld  No.l,  the  eiacoua  influence  re¬ 
mained  concentrated  only  in  the  regiona  eery  cloae  to  the 
plate  or  fln  aurfaee  (typically  upto  0.0Sda)aad  did  not  af¬ 
fect  the  (k>w  higher  up  in  the  boundary  layer.  Again,  thia 
rcault  may  be  a  conaequenee  of  inadequate  grid  reaolntioa 
for  thia  flowfleld  and  conaequently,  no  Arm  concluaiona  can 
be  drawn  regarding  the  eiacoua  influence  in  the  flow  from 
the  analyda  of  Flowfleld  No.2. 

Caae  2i  (<x,A)s(24*,40*)  Swept  Compreaaioa  Cor* 
aer  Interaction 

Aa  obaerred  from  Fig.7,  the  'rortical  atructure’  in  thia 
interaction  ia  aligned  approximately  with  the  baae  of  the 
compreadon  ramp.  In  order  to  atudy  the  interaction  and 
the  nature  of  the  rortical  atructure  ia  a  «imiler  to 

that  of  the  aharp  fln  cooe,  the  preaaure  and  eiacoua  terma 
ia  the  mean  kiaetie  energy  equation  were  compared  quaa- 
tilatirdy  for  point!  lying  on  ‘p  -  a*  planea  ia  the  flow. 
Here,  ‘y’  may  be  conddered  to  be  the  hdght  of  a  point 
aboee  the  aurfaee  of  the  plate/ramp. 

Flowfleld  No.l. 

Fig.29  ahoera  the  coatoure  for  the  iaeiadd  term  (A), 
the  eiacoua  term  (B)  and  the  preaaure  (C)  and  eelodty 
rector  (D)  plota  tor  the  *p  —  a’  planea  at  aasd.M^  and 
amT.OdeltOw.  Theae  planea  repreaent  two  'typical*  croao- 
aectiona  of  the  flow,  auflldently  (ar  from  the  origin  of  the 
rortical  atructure  aligned  approximatdy  with  the  comer. 
The  preaaure  contoura  ahow  the  diatributioa  of  the  prea¬ 
aure  rioe  dtM  to  the  ahock  orer  a  fbiile  area  aa  the  ahock 
interacta  erith  the  incoming  turbulent  boundary  layer  on 
the  flat  plate,  thua  forming  a  'X’  ahaped  region  of  preaaure 
riae,  dmilar  to  that  obaerred  tot  the  caoe  of  the  aharp  fln 
interaction.  The  region  of  maximum  inriodd  influence 


(ace  contour!  of  ‘inriadd  term’)  ia  concentrated  towarda 
the  upatream  aide  of  the  'A*  ahock  atructure.  Aa  in  the 
caae  of  the  aharp  fin  interaction,  the  ‘rcro’  inriadd  influ¬ 
ence  to  the  left  of  thia  region  maybe  due  to  the  orien¬ 
tation  of  the  rdodty  eomponeata  and  the  eorreaponding 
preaaure  gradient!  (in  the  ‘inriadd*  term)  in  thia  region. 
Alao,  a  significantly  eiacoua  region  liea  between  the  region 
of  maximum  inriadd  influence  and  the  compreadon  romp 
surface.  This  structnre  is  in  eery  dose  agreement  with 
that  obaerred  for  the  aharp  fln  configuration.  The  ‘polar 
rdodty  rector’  plot  representa  the  reaultant  rdodtica  of 
the  components  normal  to  a  ray  drawn  from  the  origin  (at 
the  baae  of  the  comer  and  at  the  asO  boundary  of  the  do- 
main).  The  directions  of  the  rdodty  rectors  shows  quite 
dearly,  the  rortical  atracturc  at  the  base  of  the  comer. 

Flowfldd  Wo.2. 

Fig.30  shows  the  contours  far  the  inriadd  and  riacous 
terms  (A,  B)  and  preaaure  (C)  along  with  the  *polar  rdoc- 
ity  plot*  (D)  for  the  'y-a*  planes  at  tas4.7t^  and  taiT.Ofm- 
The  region  of  inriadd  influence  again  lies,  to  a  large  ex¬ 
tent,  towards  the  upstream  ddc  of  the  ‘A*  ahaped  region 
of  pressure  rise.  The  region  of  riacous  influence  eon  be 
obaerred  to  lie  dong  the  atreamwiae  direction  at  a  eertdn 
distance  from  the  surface  of  the  plate  or  the  compreadon 
comer.  The  region  of  maximum  riacous  influence  liea  (ar 
downstream  of  the  region  of  concentrated  inriadd  influ¬ 
ence  in  the  flow.  Agdn,  howerer,  it  is  poedble  that  the  in> 
adequate  grid  resolution  has  introduced  dgnifleant  errors 
in  the  computation  of  the  gradients  of  the  riacous  atreaaea 
for  points  lying  on  these  planes.  Hence,  no  firm  eondu- 
dons  can  be  drarm  from  the  andyda  of  thia  flowfldd  with 
regard  to  regions  of  riacous  influence  in  the  interaction. 

The  detailed  inrestigation  of  the  riacous  and  inriadd 
influence  ia  the  3-D  turbulent  interactions  generated  by 
the  sharp  fin  and  swept  comer  conflgurations  yidds  the 
following  results: 

1.  The  dynamics  of  the  interactions  am  obaerred  to  be 
goremed  primarily  by  inriadd  (preaaum)  eflecta. 

2.  Viacoua  influence  is  dgnifleant  ia  the  interactions,  as 
obaerred  from  the  ezaminatioa  of  the  flowfldda  generated 
by  the  Narier-Stokes  code  incorporating  the  Bddwin- Lomax 
dgebrdc  turbulent  eddy  riacodty  modd.  In  the  ease 
of  the  aharp  fin,  this  region  extends  upto  approodmatdy 
Z.OSm  abore  the  plate  aurfiace. 

3.  The  region  of  maximum  inriadd  influence  ia  the  inter¬ 
actions  is  concentrated  towarda  the  freestreom  dde  of  the 
A  ahaped  region  of  preaaum  rise  (or  the  aharp  fln  inter¬ 
action  and  towards  the  upatmam  dde  for  the  swept  com- 
preadon  comer  case.  Thia  is  obaerred  to  be  condatent  for 
both  flowfldda  (generated  udng  different  turbulence  mod¬ 
els)  ia  each  case. 
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4.  The  re^oa  of  iTi««<minii  ritcoua  influence  lies  betweea 
the  region  of  maximum  inTiscid  influence  and  the  geome¬ 
try  generating  the  interaction  for  the  flosrflelds  computed 
using  the  Baldwin- Lomax  model, 

5.  Mo  Arm  eoncluaioas  can  be  drawn  regarding  region  of 
ritcoua  influence  in  the  interactions,  by  the  analytit  of  the 
flowflelda  generated  using  the  Jonea-Launder  turbulence 
model,  due  to  inadequate  grid  resolution  for  computation 
of  the  riscous  streaa  gradient  terma. 

The  present  study  gires  an  indication  of  the  complex 
nature  of  the  3-D  shock  ware-  tnrbnlent  boundary  layer 
interaction  generated  by  eren  simple  dimensionless  ge¬ 
ometries.  It  would  be  interesting  to  examine  the  flowflelda 
with  regard  to  rortidty  and  total  energy  along  streamlines 
in  the  flow.  It  would  be  also  worthwhile  to  compare  the 
results  obtained  foe  the  sharp  fln  and  swept  eomer  inter¬ 
action  to  other  three  dimensiosal  intemetiona. 
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Fig.S  Kinetic  energy  along  (treatnline 
originating  from  y  =0.024^ 

Flowfield  No.l 

Fig.9  Kinetic  energy  along  streamline 
originating  from  y  =0.05$^ 

Flowfield  No.l 

Fig.lO  Kinetic  energj  dong  itrenmline 
originating  &om  y  =0.2i. 
Floirfield  No.l 


Fig.  11  Kinetic  energy  dong  strearaline 
originating  bom  y  =0.Sd« 
Flowiield  No.l 


Fig.l2  Streamline*  for  Swept  compression  comer 
(a,A)=(24,40)  at  Res.  =  8.1x10* 

Flowfield  generated  using  the  Jones-Launder  turbulence  model 


Fig.  13  Kinetic  energy  dong  streamline 
originating  from  y.  =0.02S^ 
Flowdeld  No.2 
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Fig.l4  Kinetic  energy  dong  streamline 
originating  from  y  =0.05S^ 
Flowfield  No.2 
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Fig.15  Streamline*  for  Swept  compresiioo  cornet 
(a, A)  =(24.40)  at  Re«.  =  8.1x10* 

Floi^eld  generated  using  the  Jone*-Launder  turbulence  model 
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Fig.ld  Kinetic  energj  along  streamline 

Fig.l7  Kinetic  energy  along  streamline 

ocigittating  from  y  =0.25.. 

originating  from  j  =0.5f„ 
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Flowfield  No.2 
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Fig.lS  Streamlinea  for  tharp  nn  connguration 
(a  =20*)  at  Rei.  =  9.0x10* 

Flowfield  generated  uaing  the  Baldwin-Lomajc  turbulence  model 
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Fig.l9  Kinetic  energj  along  itreamline 

Flowfleld  No.l,  Sharp  fin  interaction 
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Fig.20  Kinetic  enetgj  along  itreamline 

Floirfielfi  No.l,  Sharp  fin  interaction 


Fig.21  Kinetic  energy  along  itreamline 

Flowfielfi  No.l,  Sharp  fin  interaction 
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Fig.23  Preiiure/Viicoiia  Contribution, Height  above  plate 
Flowfield  No.l,  Sharp  fin  interaction 
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Fig.22  Kinetic  energy  along  itreamline 

Flowfield  No.l,  Sharp  fin  interaction 
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Fig.24  Preaiure/Viicous  Contribution,Height  above  plate 
Flowfield  No.l,  Sharp  fin  interaction 
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Fig.2<  Kin<tie  energy  elong  ttreeniline 

Flowfield  No.S,  Sharp  fin  interaction 


Fig.27  Kinetic  energy  along  itreanihne 

Plowfield  No.2,  Sharp  fin  interaction 
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F:^-29  Inviidd,  dtcou*  contribution  to  rnte  of  change 
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‘Polar’  relodty  rector* 

(a,A)=  (24,40]Swept  compreiiion  comer  interaction 

Flowfield  No.l,  Baldwin- Lomax  turbulence  model 
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Fig.30  Inriacid,  riaeoua  contribution  to  rate  of  change 
of  mean  kinetic  energy,  Preaaure  contoura, 

‘Polar’  Telocity  rectora 

(a,A)=  (24,40)Swept  eompreaaion  comer  interaction 
Flowfield  No.2,  Jonea-Launder  turbulence  model 


